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ABSTRACT 

We report cosmic microwave background (CMB) power spectrum measurements from the first 
100 deg 2 field observed by the South Pole Telescope (SPT) at 150 and 220 GHz. On angular scales 
where the primary CMB anisotropy is dominant, t < 3000, the SPT power spectrum is consistent with 
the standard ACDM cosmology. On smaller scales, we sec strong evidence for a point source contribu- 
tion, consistent with a population of dusty, star-forming galaxies. After we mask bright point sources, 
anisotropy power on angular scales of 3000 < i < 9500 is detected with a signal-to-noise > 50 at 
both frequencies. We combine the 150 and 220 GHz data to remove the majority of the point source 
power, and use the point source subtracted spectrum to detect Sunyaev-Zel'dovich (SZ) power at 
2.6 cr. At t = 3000, the SZ power in the subtracted bandpowers is 4.2 ± 1.5 /iK 2 , which is significantly 
lower than the power predicted by a fiducial model using WMAP5 cosmological parameters. This 
discrepancy may suggest that contemporary galaxy cluster models overestimate the thermal pressure 
of intracluster gas. Alternatively, this result can be interpreted as evidence for lower values of og. 
When combined with an estimate of the kinetic SZ contribution, the measured SZ amplitude shifts 
erg from the primary CMB anisotropy derived constraint of 0.794 ±0.028 down to 0.773 ±0.025. The 
uncertainty in the constraint on erg from this analysis is dominated by uncertainties in the theoretical 
modeling required to predict the amplitude of the SZ power spectrum for a given set of cosmological 
parameters. 

Subject headings: cosmology - cosmology: cosmic microwave background - cosmology: observations - 
large-scale structure of universe 



1. INTRODUCTION 
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Anisotropy in the cosmic microwave background 
(CMB) has been well-char acterized on angular scales 
larger than a few arcminutes ( Jones et al. 2006 ; Rei chardt] 



et al.||2009a| INolta et al. 2009" Brown eh al. 2009b, pro 



viding stringent tests ot cosmological theory a nd strong 



evidence fo r the ACDM cosmolo gical model (Hinshaw 
et al.||2009| |Komatsu et al.|[2009[) The primary CMB 
anisotropy is a direct record of density fluctuations at 
the surface of last scattering. On angular scales < 10 
arcminutes, the primary CMB anisotropy is exponen- 
tially dampe d due to photon diffusion in the primor- 
dial plasma ( |Silk| [l968); the resulting decline in power 
with increasing multipole is known as the "damping tail" . 
The anisotropy on very small scales, which is only begin- 
ning to be explored experimentally, is instead dominated 
by foreground emission and secondary distortions to the 
CMB introduced after the surface of last scattering. This 
secondary CMB anisotropy is produced by interactions 
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of CMB photons with large scale structure and can po- 
tentially provide sensitive and independent tests of cos- 
mological parameters important to structure formation. 

The most significant source of secondary CMB temper- 
ature anisotropy on small angular scales is predicted to 
be the inverse Compton scattering of CMB photons by 
the hot plasma gravitationally bound to m assive, col- 
lapsed objects (Sunyaev & Zeldovich 1972[ ). This is 
known as the thermal Sunyaev-Zel'dovich (tSZ) effect. 
The tSZ effect distorts the blackbody spectrum of the 
primary CMB, creating a decrement in intensity at low 
frequencies and an increment at high frequencies with 
a null near 220 GHz. At 150 GHz, the anisotropy in- 
duced by the tSZ effect is expected to dominate over 
the primary CMB fluctuations at multipoles £ > 3000. 
The power contributed to CMB maps by the tSZ effect 
depends sensitively on the normalization of the matter 
power spectrum, as parametrized by the RMS of the mass 
distribution on 8h _1 Mpc scales, <j%. 

Galaxy clusters over a wide range of mass and redshift 
contribute significantly to the tSZ power spectrum. It 
is challenging to accurately model the low-mass or high- 
redshift clusters. For instance, lower mass clusters are 
more sensitive to non-gravitational heating effects such 
as by AGN, while we have limited observational data 
on high-redshift clusters. There is considerable uncer- 
tainty as to the expected shape and amplitude of the 
tSZ power spectrum due to the current lack of knowl- 
edge of the properties of intracluster gas in low mass and 
high redshift galaxy clusters. 

The scattered CMB photons also obtain a net Doppler 
shift when ionized matter is moving with respect to the 
rest frame of the CMB. This so-called kinetic Suny aev- 
Zel'dovich (kSZ) effect ( [Sunyaev fe Zeldovich||l980| de- 
pends solely on the motion and density of free electrons. 
In contrast to the tSZ effect, the kSZ effect has contribu- 
tions from electrons with temperatures as low as 10 4 K. 
Therefore higher-redshift epochs, before massive objects 
finish collapsing, are expected to have relatively larger 
contributions to the kSZ power. Recent simulations and 
analytic models also predict a sizable signal from the 
epoch of the first radiative sources which form ionized 
regions several tens of Mpc across, within a largely neu- 
tral Universe. Low-redshift galaxy clusters dominate the 
power on small angular scales, while high-redshift reion- 
izing regions have their largest relative contribution on 
angular scales around £ — 2000. At 150 GHz, the kinetic 
effect is expected to amount to tens of percent of the 
total SZ power. 

Only a handful of experiments have had sufficient sen- 
sitivity and angular resolution to probe the damping tail 
of the C MB anisotropy. Early measurements a t 30 GHz 



sky means that there is no significant tension between 
the BIMA and SZA results. The latest CBI measure- 



by CBI (Mason et al. 



2003| |Bond et al.]|2005| ) reported 



a > 3 a excess above the expected CMB power at multi 
poles of I > 2000. Observatio ns with the BIMA array at 
30 GHz ( Dawson et al. 2006[ ) also reported a nearly 2 a 
detection of excess power at £ = 5237. However, more 
recently, the SZA experiment (also observing at 30 GHz) 
has published an upper limit of 149 /iK 2 at 95% confi- 



dence on excess power at these multipoles (Sharp et al. 
2009 1 in apparent conflict with the previous CBI and 



ments (Sievers et al. 



2009 1 include more data, improved 



radio source removal, and a proper treatment of non- 
Gaussianity of the SZ sky. These measurements continue 
to suggest excess power but with a significance of only 

1.6a. 

At 150 GH z, the ACBAR ( iReichardt et al.||2009a l 
and QUaD ( Friedman et al.| |2009| experiments have 
both measured the damping tail of the primary CMB 
anisotropy at £ < 3000 with high signal to noise. Either 
with or without the addition of the expected foreground 
and tSZ contributions, the power measured at the highest 
multipoles by both experiments is consistent with pri- 
mary CMB anisotropy. In the last year, the results of 
150 GHz observations out to t = 10000 made with the 



Bolocam (Bayers et al.||2009| ) and APEX-SZ (Reichardt 
et al.||2009bp experiments have been released. These ex- 
periments have been used to place upper limits on power 
above the primary CMB of 1080 fjK 2 and 105 ^K 2 re- 
spectively at 95% confidence. The constraints on erg from 
these upper limits remain weak, in no small part due to 
the large, highly non-Gaussian sample variance of the 
tSZ effect on the small ~1 deg 2 patches of sky observed 
by Bolocam and APEX-SZ. The cosmic variance of the 
tSZ effect will be significantly reduced in the on-going 
> 100 deg 2 surveys being conducted by next-generation 
experiments such as ACT (?) and SPT. 

In addition to the tSZ and kSZ effects, foreground 
emission is important on these small angular scales. Af- 
ter bright radio sources are removed, the most significant 
foreground at 150 and 220 GHz is expected to be a pop- 
ulation of unresolved, faint, dusty, star-forming galax- 
ies (DSFGs) with a rest frame emission spectrum that 
peaks in the far infrared. These sources have been stud- 
ied at higher frequen cies close to the peak of their emis 



sion spectrum 21 (e.g 



Holland et al. (1999); Kreysa et al 



( [1998| ; |Glenn et al.| ( |1998| ); |Viero et al.| ( |2009| ), however 
extrapolating their fluxes to 150 GHz remains uncertain 
Adding to the challeng e is the expected significant clus- 
tering of these sources dH aiman & Knox 2000 ; Knox et al. 
jRighi et al. |2008j [Sehgal et al.|f2TO| ). IU eniis- 



2001 ^ 

sion from clustere 



DSFGs was first observed with the 
Spitzer telescope at 160 (im ( [Lagache et al. 2007) and 



BIMA results. For the relatively small patch (0.1 deg ) 
observed by BIMA, the non-Gaussian nature of the SZ 



more recently this clustering has also been obser ved at 
sub-mm wa velengths by the BLAST experiment (Viero 
et al. |[20Q9] ) . The clustering of these DSFGs is expected 
to produce anisotropic power at 150 GHz with an an- 
gular power spectrum that is similar to that of the SZ 
effect. However, emission from DSFGs is spectrally sep- 
arable from the SZ effect and the SZ power spectrum can 
be recovered by combining information from overlapping 
maps at 150 and 220 GHz. 

In this work, we present measurements by the South 
Pole Telescope (SPT) which comprise the first significant 
detections of anisotropy power for £ > 3000 at 150 and 
220 GHz. The SPT has sufficient angular resolution, sen- 
sitivity and sky coverage to produce high-precision mea- 
surements of anisotropy over a range of multipoles from 
~ 100 < £ < 9500. However, for the immediate goal 

21 Since these sources are typically brightest at sub-millimeter 
wavelengths they are also referred to in the literature as sub- 
millimeter galaxies (SMGs). 
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of measuring secondary CMB anisotropics, we start with 
the first bandpower at £ — 2000 where primary CMB still 
dominates the power spectrum. We combine bandpowers 
from two frequencies to minimize the DSFG contribution 
and produce the first significant detection of the SZ con- 
tribut ion to the CMB p ower spectrum. In a companion 
paper 



Hall et al. (2009 



^ hereafter H09), the SPT power 
spectra are used to place constraints on the amplitude of 
the Poisson and clustered components of DSFG power. 
This is the first detection of the clustered DSFG power 
at mm- wavelengths. We use three results from H09 in 
setting limits on the residual po int source power in the 
DSFG-subtracted bandpowers (j |6.2| : an argument that 
residual clustered components are negligible, a measure- 
ment of the Poisson point source power at 150 GHz, and 
a limit on the scatter in spectral index between DSFGs. 

We describe the instrument, observations, beams, and 
calibration strategy in f|2] The time-ordered data (TOD) 
filtering and map-making algorithm is outlined in Sj3] 
along with the procedure to derive bandpowers from 
maps. The results of tests for systematic errors applied 
to the SPT data are discussed in ffi] and the expected 
astronomical foregrounds are described in ^j5] The band- 
powers and cosmological interpretation are given in Sj6] 

2. INSTRUMENT AND OBSERVATIONS 

The SPT is an off-axis Gregorian telescope with a 10-m 
diameter primary mirror located at the South Pole. The 
telescope is optimized to perform high resolution surveys 
of low surface brightness sources. The first receiver to 
be placed on the SPT has the primary goal of identify- 
ing a mass-limited sam ple of galaxy clusters , the fi rst of 

, here- 
cmcnt 



which were reported in Staniszewski et al. (2009 
after S09). The receiver is equipped with a 9b0-e 
array of superconducting transition edge sensor bolome- 
ters. The detectors are split between three frequency 
bands centered at 95 GHz, 150 GHz, and 220 GHz, al- 
lowing the separation of the tSZ effect from the primary 
CMB anisotropy and foregrounds. Fu rther details of th e 
receiver and telesco pe c an be found inFRuhl et aL ( 2004 1 , 
|Padin et al.] (|2008[ ) and |Carlstrom et aTTf2TjMj K 

In this work, we use data at 150 and 220 GHz from one 
100 deg 2 field centered at right ascension 5 h 30 m , decli- 
nation -55° (J2000) observed by SPT in the first half 
of the 2008 austral winter. Data from SPT's third (95 
GHz) frequency band was unusable in 2008. The loca- 
tion of the field was chosen for overlap with the Blanco 
Cosmology Survey (BCS) 22 optical survey and low dust 
emission. We observed this field for a total of 779 hours, 
of which 575 hours is used in the analysis after passing 
data quality cuts. The final map noise is 18 ^K-arcmin 23 
at 150 GHz and 40 /JK-arcmin at 220 GHz. This field ac- 
counts for half the sky area observed in 2008 and one 
eighth of the total area observed by SPT to date. 

The scanning strategy used for these observations in- 
volves constant-elevation scans across the 10° wide field. 
After each scan back and forth across the field, the tele- 
scope executes a 0.125° step in elevation. A complete 

22 http: / /cosmology.illinois.edu/BCS 

23 Throughout this work, the unit K refers to equivalent fluc- 
tuations in the CMB temperature, i.e., the level of temperature 
fluctuation of a 2.73 K blackbody that would be required to pro- 
duce the same power fluctuation. The conversion factor is given by 
the derivative of the blackbody spectrum, evaluated at 2.73 K. 



set of scans covering the entire field takes approximately 
two hours, and we refer to each complete set as an obser- 
vation. Successive observations use a series of different 
starting elevations to ensure even coverage of the field. 
Of the 300 observations used in this analysis, half were 
performed at an azimuth scanning speed of 0.44° per sec- 
ond and half at a speed of 0.48° per second. 

2.1. Beam Measurements 

The SPT beams are measured by combining maps of 
three sources: Jupiter, Venus, and the brightest point 
source in the 100 deg 2 field. Observations of Jupiter are 
used to measure a diffuse, low- level sidelobe in the range 
15' < r < 40', where r is the radius to the beam center. 
Although this sidelobe has a low amplitude (— 50 dB at 
r=30') , it contains approximately 15% of the total beam 
solid angle. A measurement of this sidelobe is neces sary 



for the cross-calibration with WMAP described in S2.2 



The observations of Jupiter show signs of potential non- 
linearity in the response of the detectors for r < 10'. For 
this reason we only use the observations of Jupiter to 
map the sidelobe at r > 15'. Observations of Venus are 
used to measure the beam in the region 4' < r < 15'. 
The angular extent of Jupiter or Venus has a negligible 
effect on the measurement of the relatively smooth beam 
features present at these large radii. The brightest point 
source in the map of the 100 deg 2 field is used to measure 
the beam within a radius of 4'. In this way, the random 
error in the pointing reconstruction (7" RMS) and its 
impact on the effective beam are taken into account. The 
pointing error has a negligible effect on the relatively 
smooth outer (r > 4') region of the beam. 

A composite beam map, B{9 1 </>), is assembled by merg- 
ing maps of Jupiter, Venus, and the bright point source. 
The TOD are filtered prior to making these maps, in 
order to remove large-scale atmospheric noise. Masks 
with radii of 40', 25', and 5' are applied around the lo- 
cations of Jupiter, Venus, and the point source, respec- 
tively. These masks ensure that the beam measurements 
are not affected by the filtering. 

Using the the flat-sky approximation, we calculate the 
Fourier transform of the composite beam map, B(£, fy). 
From this, we compute the azimuthally-averaged beam 
function, 



B t = Re{ 



1 

2^ 



B(e, 0,)dM. 



(i) 



We note that averaging \B(£,(f>£)\ 2 instead of B(l,<j>i) 
would result in a percent-level noise bias in Be at very 
high multipoles due to the presence of noise. The results 
in this work assume an axially symmetric beam, which is 
only an approximation for SPT. We simulate the effects 
of ignoring the asymmetry on the bandpowers and find 
that the errors introduced by making this assumption are 
negligible. 

Although the measured beam function Bg is used for 
the bandpower estimation, an empirical fit is used to 
quantify the errors on Bg. Bi is fit to the empirical model 

B e = ae-i^ 5 + (1 - a ) e -i(o.oaxa«) 1 -\ (2) 

There are two components: a main lobe (first term) and 
a diffuse shelf (second term). The form of the model 
and the numerical values of the slopes of the exponents 



4 




0.10 



0.08 



0.06 



0.04 



- 0.02 



0.00 



2000 



4000 6000 



8000 10000 



Fig. 1. — Average beam functions and uncertainties for SPT. 
Left axis: The SPT beam function for 150 GHz (red) and 220 GHz 
(blue). Right axis: The lrr uncertainties on the beam function 
for each frequency. The beam uncertainties shown here include 
only uncertainties on the main lobe beam width, tr;,, since the 
uncertainty of the sidelobe amplitude has been subsumed into the 
calibration uncertainty. 

were constructed to provide a good fit to the measured 
Bi . We note that Bj 50 and B^ are measured and fitted 
separately. The RMS difference between the model and 
measured Bgs is approximately 1%. Two parameters re- 
main free: ab, which describes the width of the main lobe, 
and a, which sets the relative normalization between the 
main lobe and the diffuse shelf. These parameters are left 
free to quantify the uncertainty in Bg. The uncertainty 
in the values of these parameters directly translates to 
an uncertainty in Bg. 

There are a number of factors that limit the accu- 
racy of the measurement of Bg. These include residual 
map noise, errors associated with the map-merging pro- 
cess, and spectral differences between the CMB and the 
sources used to measure the beam. The final uncertain- 
ties in the beam model parameters a b and a are con- 
structed as the quadrature sum of the estimated uncer- 
tainties due to each of these individual sources of error. 
To a good approximation, the uncertainties on <7b and a 
can be taken to be Gaussian and uncorrelated. 

In practice, a change in the value of a is equivalent 
to a change in the overall calibration for I > 700. Af- 
ter the beam uncertainties are estimated, the uncertainty 
in a is folded into the estimated uncertainty on the ab- 
solute calibration and the parameter a is fixed to the 
best-fit value of 0.8 5. The quoted beam uncertainties 
in the CosmoMC 24 ( |Lewis fc Bridle||2002[ ) data files in 
|6.1| include only the uncertainty on Ob- Figure [T] shows 
the measured beam functions for 150 and 220 GHz, along 
with the 1 a uncertainties in the main lobe beam width 
a b . 

2.2. Calibration 

The calibration of the SPT data is tied to the superb 
WMAP5 absolute calibration through a direct compari- 
son of 150 GHz SPT maps with WMAP5 V and W-band 

24 http:/ /cosmologist.info/cosmomc 



(61 and 94 GHz) maps ffiinshaw et al.|2009[ ) of the same 
sky regions. The WMAP5 maps are resarripled accord- 
ing to the SPT pointing information, and the resulting 
TOD are passed through the SPT analysis pipeline to 
capture the effects of TOD filtering. The ratio of the 
cross-spectra of the filtered WMAP and SPT maps after 
correcting for the instrumental beams, 



a lm,WMAPi a tm,WMAPj 



l lm,SPT a £m,WMAP j 



(3) 



is used to estimate the relative calibration factor between 
the two experiments. A similar procedure was used to 
calibra te the Boomerang, ACBAR, and QUaD experi- 
ments Jjjones et ai]|2006| Reichardt et al.||2009a| |Brown" 
et al. 2009p . Dedicated SPT calibration scans of lour 
large fields totaling 1250 deg 2 of sky were obtained dur- 
ing 2008. The results for these four fields are combined to 
achieve an absolute temperature calibration uncertainty 
of 3.6% at 150 GHz. 

The 150 GHz calibration is transferred to 220 GHz 
through the overlapping coverage of SPT's high S/N 
maps. We calculate the relative calibration by exam- 
ining the ratio of the cross-spectra between the 150 and 
220 GHz maps to the auto-spectra of the 150 GHz map 
after correcting for the beam and filtering differences. 
We estimate the relative calibration uncertainty to be 
6.2% and the final absolute calibration uncertainty of 
the 220 GHz temperature map to be 7.2%. 

3. ANALYSIS 

In this section we describe the pipeline used to process 
the TOD to maps, and from maps to bandpowers. Given 
the small sky area analyzed here, we use the flat-sky ap- 
proximation. We generate maps in a flat-sky projection 
and analyze these maps using Fourier transforms. Thus, 
the discussion of filtering and data-processing techniques 
refer to particular modes by their corresponding angu- 
lar wavenumber k in radians. Note that in these units 
\k\=£. 

3.1. Map-making 

We generate one map per frequency band for each two- 
hour observation of the field. The map-making pipeline 
used for this analysis is very similar to that used by 
S09, with some small modifications. An overview of the 
pipeline is provided below, emphasizing the aspects of 
data processing and map-making that are most impor- 
tant for understanding the power spectrum analysis. 

3.1.1. Data Selection 

Variations in observing conditions and daily receiver 
setup can affect the performance of individual bolometer 
channels. The first step in the data processing is to select 
the set of bolometers with good performance for each in- 
dividual observation. The criteria for selecting these de- 
tectors are primarily based on responsivity (determined 
through the series of calibrations performed prior to each 
field observation) and noise performance. Detectors are 
rejected if they have a low signal-to-noise response to a 
chopped thermal calibration source or atmospheric emis- 
sion (modulated through a short elevation scan). They 
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are also rejected if their noise power spectrum is heavily 
contaminated by readout-related line-like features. More 
detail on the bolometer cuts can be found in S09. In 
addition, bolometers that have a responsivity or inverse- 
noise-based weight more than a factor of three above or 
below the median for their observing band are rejected. 
On average, 286 bolometers at 150 GHz (out of 394 possi- 
ble) passed all cuts for a given observation. At 220 GHz, 
an average of 161 (out of 254 possible bolometers) passed 
all cuts. 

In a given observation, segments of the TOD corre- 
sponding to individual scans may be rejected due to 
readout or cosmic-ray induced features, anomalous noise, 
or problems with pointing reconstruction. The scan-by- 
scan data cuts are the same as those described in S09 
and remove approximately 10% of the TOD. 

The SPT detectors exhibit some sensitivity to the re- 
ceiver's pulse-tube cooler. Bolometer noise power spectra 
are occasionally contaminated by narrow spectral lines 
corresponding to the pulse-tube frequency and its har- 
monics. As in S09, we address this by applying a notch 
filter to the affected data. This filter removes less than 
0.4% of the total signal bandwidth. 

In addition to the cuts above, we remove a small num- 
ber of observations that have either incomplete coverage 
of the field or high noise levels. Observations taken un- 
der poor atmospheric conditions are rejected in this way. 
The application of these cuts remove 36 out of the 336 
total observations of the field. 

3.1.2. Time- Ordered Data (TOD) Filtering 

Let d a j be a measurement of the CMB brightness tem- 
perature by detector j at time a. Contributions to d a j 
are the celestial sky temperature, s a j, the atmospheric 
temperature, a a j, and instrumental noise, n a j. The in- 
strumental noise is largely uncorrelated between detec- 
tors. However, the atmospheric signal is highly corre- 
lated across the focal plane due to the overlap of individ- 
ual detector beams as they pass through turbulent layers 
in the atmosphere. 

The signals in the data have been low-pass filtered by 
the bolometers' optical time constants. These single-pole 
bolometer transfer functions are measured by the method 
described in S09. We deconvolve the transfer functions 
for each bolometer on a by-scan basis, as the first step 
in TOD filtering. At the same time, we apply a low- 
pass filter with a cutoff at 12.5 Hz to remove noise above 
the Nyquist frequency associated with the final map pix- 
elation. These operations can be represented as linear 
operation II fft on the TOD, 



(4) 



where a sum is implied over repeated indices. 

To remove 1 /f noise and atmospheric noise in the scan 
direction, we project out a 19th order Legendre polyno- 
mial from the TOD for each detector on a scan-by-scan 
basis. This operation 11* can be described by 



A" — TT* A 1 



(5) 



The effect on the signal is similar to having applied a 
one-dimensional high-pass filter in the scan direction at 
£ > 300. 



Significant atmospheric signal remains in the data af- 
ter this temporal polynomial removal. Because the at- 
mospheric power is primarily common across the entire 
focal plane, we can exploit the spatial correlations to re- 
move atmosphere without removing fine-scale astronom- 
ical signal. The method used in this analysis is to fit for 
and subtract a plane, aQ + a x x. + a y y, across all detectors 
in the detector array at each time sample, where x and 
y are set by the angular separation of each pixel from 
the boresight of the telescope. All detectors are equally 
weighted in this fit. This operation can be represented 
as a linear operation on all detectors at each time sample 
a to produce an atmosphere-cleaned dataset, d'", 

<Q = n5A- (6) 

We mask the brightest point sources in the map be- 
fore applying the polynomial subtraction and the spatial- 
mode filtering. The 92 sources that were detected above 
5 a in a preliminary 150 GHz map have been masked. A 
complete discussion of the point sources detected in this 
field can be found in jVieira et al.| ( |2009[ hereafter V09) . 

3.1.3. Map-making 

The data from each bolometer is inverse noise weighted 
according to the calibrated, pre-filtering detector PSD in 
the range 1-3 Hz, corresponding to 1400 < I < 4300. 
This multipolc range covers the angular scales where we 
expect the most significant detection of the SZ effect. We 
represent the mapping between time-ordered bolometer 
samples and celestial positions with the pointing matrix 
L, which we apply to the cleaned TOD to produce a map 



LA w n s n*n fft d 



(7) 

A w is a diagonal matrix encapsulating the detector 
weights. Informatio n on the pointing recons truction can 
be found in S09 and iCarlstrom et al.l d2009 1. 



3.2. Maps to Bandpowers 

We use a pseudo-CV method to estimate the bandpow- 
ers. In pseudo-CV methods, bandpowers are estimated 
directly from the Fourier transform of the map after cor- 
recting for effects such as TOD filtering, beams, and fi- 
nite sky coverage. We p rocess the data using a cross 



spectrum b ased analysis (Polenta et al. 2005 Tristram 



et al. 2005) in order to eliminate noise bias. Beam and 



filtering effects are corrected for according to the fo rmal- 
ism in the MASTER algorithm ( |Hivon et al~||2002| ). We 
report the bandpowers in terms of 'Vi, where 



2tt 



(8) 



As the first step in our analysis, we calculate the 
Fourier transform, mt, of the single-observation maps 
for each frequency. Each map is apodized using the same 
window W, and thus fh^' A ^ = FT (Wm' 1 ''' 4 )) , where the 
first superscript, v, indicates the observing frequency, the 
second superscript, A, indicates the observation number, 
and the subscript denotes the angular frequency. Cross- 
spectra are computed for every map pair, and averaged 
within the appropriate (,-hin b: 



D 



-E 



k(k+ 1) _ f u A \ _ ( Vj ,B)* 



2tt 



(9) 



We use the abbreviation, D Ui = D UiXUi , when referring 
to single- frequency auto-spectra. Recall that |k| = £ in 
the flat-sky approximation. 

The FFT of a map, m^' , is linearly related to the 
FFT of the astronomical sky, a£. It also includes a noise 



contribution, 



(v,Ai) 



which is uncorrelated between the 



300 observations: 



B 



k'« k ' 



(10) 



Here Gk is the 2-dimcnsional amplitude transfer func- 
tion, which accounts for the TOD fi ltering as well as 

the map-based filtering described in §3.2.2 is the 



Fourier transform of the apo dization mas k 



Following the treatment by Hivon et al. (2002 1, we take 



the raw spectrum, D^ B , to be linearly related to the true 
spectrum by a transformation, . The K transforma- 
tion combines the power spectrum transfer function, F^, 
which incl udes t he effects of TOD- base d and map-based 
filtering ({ 3.2.3), the beams, Bi 
coupling matrix, M kk >[W] ( §3.2.1 1 



2.1 



and the modc- 
le mode-coupling 



matrix accounts for the convolution of the spectrum due 
to the apodization window, W. 



Kv"' = p bk ( M kk ,[W] F v k j^B%B v j ) Q Vh ,. (11) 



Here, Pb k is the re-binning operator ( Hivon et al.||2002 l 



PI 



bk 



fc(fc+l) 
2nAk 




(i>) 



fe(6_i) < k < k 



0>) 



otherwise 



while the inverse of the re-binning operator is Q 



kb- 



!kb 



fc(fri) hb-i) < k < k m 

otherwise 



(12) 



(13) 



In the rest of the paper, we will often refer to band aver- 
aged quantities, such as C& = Pb k C k . Cross-spectra be- 
tween observations that have been corrected for apodiza- 
tion and processing are denoted as, 

The final bandpowers are then computed as the average 
of all cross-spectra, 



,AB 



n-obs (n-obs - 1) 



(15) 



A B^A 



We make the simplifying assumption that the noise 
properties of each observation are statistically equiva- 
lent, hence the uniform weigh ting chosen here. The 



data selection criteria in £3.1.1 ensure that these obser- 



vations target the same region of the sky, have roughly 
the same number of active bolometers, and have simi- 
lar noise properties. Therefore, this uniform weighting 
should be unbiased and only slightly sub-optimal. 

The cross-spectrum bandpowers, D^' x Vi ' B , generated 
from the 300 observations are used in conjunction with 
signal-only Monte Carlo bandpowers to generate empiri- 
cal covariance matrices, as described in appendix [X] The 
variance of the Monte Carlo bandpowers is used to esti- 
mate the sample variance contribution to the covariance 
matrix. Meanwhile, we use the variance in the spectra 



of the real data to estimate the uncertainty due to noise 
in the maps. 

3.2.1. Apodization Mask and Calculation of the 
Mode-mixing Kernel 

Since we have only mapped a fraction of the full sky, 
the angular power spectra of the maps are convolutions 
of the true C^s with an ^-space, mode-mixing kernel 
that depends on the map size, apodization, and point 
source masking. We calculate this mode-mixing kernel , 
M kk i [W] , following the derivation in Hivon et al. ( 2002 1 
for the flat sky case: 



M kk ,[W] 



1 



(27T) 



k-k' 



(16) 



If the mask is smooth on fine angular scales, then the 
mode coupling kernel can be approximated by a delta 
function at high-fc: 

M kk ,[W]^w 2 S, 



kk>- (17) 

Here we use the notation w n = (W") to represent the 
n th moment of the apodization window. In this limit, 
the coupling kernel serves the purpose of re-normalizing 
the power spectrum to account for modes lost due to 
apodization. As we discuss in appendix [Aj the coupling 
kernel also plays an important role in determining the 
shape and normalization of the covariance matrix. 

We avoid areas of the map with sparse or uneven cov- 
erage in any single observation. Thus, the apodization 
window is conservative in its avoidance of the map edges. 
We also mask 144 point sources detected in the 150 GHz 
data above a significance of 5 a, which corresponds to 
a flux of 6.4 mjy This source list is from a more re- 
fined analysis than the preliminary list used in §3.1. 2| the 
differences are in sources near the 5 a detection thrcsh- 
old. Each point source is masked by a 2'-radius disk, 
with a Gaussian taper outside this radius. Many differ- 
ent tapered mask shapes were tested for both efficacy 
in removing point source power and noise performance. 
Given the relatively small area masked by point sources, 
varying the shape of the point source mask has little ef- 
fect on the final spectrum. The effective sky area of the 
mask is78 deg 2 . Simulations have been performed to test 
whether the application of this mask will bias the in- 
ferred power and we see no bias with a 5 a cut. As could 
be expected, we do observe a mild noise bias when using 
a more aggressive 3 a point source mask. The same 5 a 
mask is used for all maps at both frequencies and for all 
observations. 

3.2.2. Fourier Mode Weighting 



The maps produced by the steps described in §3.1.2| 
have anisotropic signal and noise. In particular, the map 
noise PSD rises steeply at spatial frequencies correspond- 
ing to low spatial frequencies in the scan direction (low 
k x ). The covariance of the power estimated in a given 
£-hva depends on the second power of the map noise PSD 
for all the Fourier modes in that £-bin. In the presence of 
either non-uniform noise or signal, applying an optimal 
mode-weighting when calculating the mean bandpower 
may significantly reduce the final noise covariance matrix 
of the power spectrum bandpowers. In the case of SPT, 
we find that for the purposes of measuring the £ > 2000 
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power spectrum, a simple, uniform selection of modes at 
k x > 1200 is close to optimal, and we apply this mode 
weighting in calculating the SPT bandpowers. 

3.2.3. Transfer Function Estimation 

In order to empirically determine the effect of both 
the TOD-based filtering and the Fourier mode- wei ghting, 
we cal culate a transfer function, F k , as defined in |Hivon| 
( 2002 1 . Note that this power spectrum transfer 



initial guess of the transfer function: 



et al 



function is distinct from the amplitude transfer function 
Gk- Specifically, the transfer function accounts for all 
map- and TOD-processing not taken into account by the 
mode-coupling kernel or the beam functions. 

We created 300 Monte Carlo sky realizations at 150 
and 220 GHz in order to calculate the transfer function 
of the filtering. The simulations also serve as an input 
for the covariance matrix estimation. These simulations 
contain two components: a CMB component and a point 
source component. The CMB component is computed 
for the best-fit WMAP5 lensed ACDM model. 

The point source component includes two different 
populations of dusty galaxies, a low-z population and 
a high-z population. For each population we generate 
sources from a Poisson distribution. Sources are gener- 
ated in bins of flux, S, ranging from 0.01 to 6.75 mJy 
This upper limit in flux is close to the 5 a detection 
threshold in the 150 GHz maps. In each flux bin, the 
150 GHz source density dN/dS, of each po pulation is 
taken from the model of Negrello et al. ( 2007 1 . We relate 
the flux of each source at 220 GHz to its flux at 150 GHz 
with a power law in intensity, S v cx v a . The power law 
spectral index, a, for each source is drawn from a normal 
distribution. We use spectral indices of a — 3 ± 0.5 for 
the high-z protospheroidal galaxies, and a = 2 ± 0.3 for 
the low-z IRAS-like galaxies. As with any Poisson dis- 
tribution of point sources, the power spectrum of these 
maps is white (C^ s = constant) and related to the flux 
cutoff, 5*0, by: 



cr 



dBu 
dT 



Tern 



So dN 

S 2 ^dS. 
dS 



(18) 



The power spectra of these simulated point source maps 
are well represented by a constant Cf s = 1.1 x 10~ 5 /iK 2 
at 150 GHz and Cf = 6.8 x 10~ 5 ^K 2 at 220 GHz. 

These simulated maps are smoothed by the appropri- 
ate beam. From each map realization, we construct sim- 
ulated TOD using the pointing information. Each real- 
ization of the TOD is processed using the low-pass filter, 
polynomial r emova l, and the spatial-mode subtraction 
described in §3.1.2| No time-constant deconvolution is 
applied, since these realizations of the TOD have not 
been convolved by the bolometer time constants. The 
filtered, simulated TOD are then converted into maps 
according to equation [7] 

We apply the same apodization mask and Fourier mode 
weighting to these map realizations as is used on the 
actual data. We then compute the Monte Carlo pseudo- 
power spectrum, (Dfc) MC , for each map. The simulated 
transfer function is calculated iteratively, by comparing 
the Monte Carlo average, (Dk) MC , to the input theory 

i ^theory 

spectrum, C k 

For the single-frequency bandpowers, we start with an 



F, 



",(0) 



w 2 Bfcf hcmy ' 



(19) 



The superscript, (0), indicates that this is the first iter- 
ation in the transfer function estimates. For this initial 
guess, we approximate the coupling kernel as largely di- 
agonal as in equation 1 17| Thus, the factor W2 is the nor- 
malization factor required by the apodization window. 
We then iterate on this estimate using the full mode- 
coupling kernel: 



= F, 



MC 



M kk ,F^B^Cl) 



2 theory 
W 



BfC^ thcmy w 2 



(20) 

We iterate on this estimate five times, although the trans- 
fer function has largely converged after the first iteration. 

This method may misestimate the transfer function if 
the simulated spectrum is significantly different from the 
true power spectrum. The primary CMB anisotropy has 
been adequately constrained by previous experiments, 
however, the foreground power spectrum is less well 
known. We repeated the transfer function estimation 
using an input power spectrum with twice the nominal 
point source power. The resulting transfer function was 
unchanged at the 1% level, giving confidence that the 
transfer function estimate is robust. 

The transfer function for a multifrequency cross- 
spectrum is taken to be the geometric mean of the two 
individual transfer functions: 



F 



F k F k 



(21) 



This treatment is only strictly correct for isotropic fil- 
tering. As a cross check, we have also computed the 
cross-spectrum transfer function directly. For the angu- 
lar multipoles reported in this work, the geometric-mean 
transfer function estimate is in excellent agreement with 



the estimate obtained using (D^ 1 
andiCD 



) MC in equations 



19 



3.2.4. Frequency-differenced Spectra 

We are interested in the power spectra of linear com- 
binations of the 150 and 220 GHz maps designed to re- 
move astronomical foregrounds. One method for obtain- 
ing such power spectra would be to directly subtract the 
maps after correcting the maps for differences in beams 
or processing. In such a map subtraction, the scaling of 
the 220 GHz map, x, can be adjusted to optimally re- 
move foregrounds. The differenced maps can then be 
processed using our standard pipeline to obtain spectra 
with a reduced foreground contribution. 

Equivalently, the differenced spectrum can be gener- 
ated from the original bandpowers, q l bl using a linear 
spectrum transformation, £: 



(22) 



Here, the index, i, denotes the 150 GHz auto-spectrum, 
220 GHz auto-spectrum, or 150 x 220 GHz cross spec- 
trum. This transformation is computationally fast and 
takes advantage of the fact that the cross-frequency 
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bandpowers include information on the relative phases 
of each Fourier component in the map. In this way, the 
difference spectrum can be represented in terms of the 
three measured spectra: 



lb - ~Z 72 l°k _ ^"k 



v 3 |2 



(1 — x) v ' 



(23) 



The overall normalization is chosen such that the CMB 
power is unchanged in the subtracted spectrum. For clar- 
ity, we have momentarily avoided here the complications 
of beams and filtering and have expressed the bandpow- 
ers in terms of the Fourier transform of the celestial sky, 
at- For a given proportionality constant, x, the values 
of £(x) are: 



-150x220 



(X): 



(l-x) 2 ' 

-2x 
(l-x) 2 ' 
x 2 

(l-x) 2 ' 



(24) 



We also compute the covariance matrix for the sub- 
tracted spectrum from the original bandpower covariance 
matrices: 



^bb' 



150-xx220 



= j2e(x)c^Pe(x). 



4. JACK-KNIVES 

We apply a set of jack-knife tests to the SPT data to 
search for possible systematic errors. In a jack-knife test, 
the data set is divided into two halves, based on features 
of the data associated with potential sources of system- 
atic error. The two halves are differenced to remove any 
astronomical signal, and the resulting power spectrum 
is compared to zero. Significant deviations from zero 
would indicate either a systematic problem or a noise 
misestimate. We implement the jack-knives in the cross- 
spectrum framework by differencing single pairs of obser- 
vations and applying the cross-spectrum estimator out- 
lined in §3.2| to the set of differenced pairs. In total, we 
perform 13 different jack-knife tests. 

Six jack-knives are based on the observing parameters, 
such as time, scan direction and azimuthal range. The 
data can be split based on when it was taken to search 
for systematic changes in the calibration, beams, detec- 
tor time constants, or any other time variable effect. The 
"first half - second half jack-knife probes variations on 
month time scales, while an "even -odd" jack-knife dif- 
ferencing every other observation looks for variations on 
hourly time scales. Results for the "first half - second 
half jack-knife are shown in the top panel of Figure [2j 
The data can also be split based on the direction of the 
scan in a "left -right" jack-knife (panel 2 of Figure [2]) . 
We would expect to see residual power here if the detec- 
tor transfer function had been improperly de-convolved, 



if the telescope acceleration at turn-arounds induces a 
signal through sky modulation or microphonics, or if the 
wind direction is important. We observed this field at 
two scan speeds different by 10%. We check for sys- 
tematic differences related to the scan speed, such as a 
mirror wobble, in a "scan speed" jack-knife. Side-lobe 
pickup could potentially introduce spurious signals into 
the SPT maps from the moon or features on the ground. 
We test for moon pickup by splitting the data based on 
whether the moon is above or below the horizon. We test 
for ground pickup by splitting on the mean azimuth of 
the observation. To maximize the sensitivity to ground 
pickup, the azimuthal ranges are selected to be centered 
on and directly opposite the closest building to the tele- 
scope, which is the most likely source of ground signal. 
The azimuthal jack-knife is shown in the third panel of 
Figured 

We also perform jack-knives based on four noise and 
observation-quality statistics of the 150 GHz data. The 
first is based on the overall RMS in the maps, which is 
affected by atmospheric conditions and detector noise. 
The second is based on the average raw detector PSD 
in the range 9-11 Hz, which is a measure of the detector 
"white" noise level. The third is the RMS near I = 3000 
where the S/N on the SZ power spectrum is highest. 
The fourth is based on the number of bolometers active 
in each observation. 

There are also a number of line-like spectral features 
in the SPT TOD that could potentially affect the power 
spectrum bandpowers, and we perform three jacknife 
tests for sensitivity to these features. Some of these 
line features appear at harmonics of the receiver pulse- 
tube frequency, and are typically correlated across many 
bolometer channels. These lines are filtered from the 
data as described in §3.1.2| In addition, some channels 
exhibit occasional line-like features at other frequencies, 
which are not filtered in the data processing. We search 
for residual effects in a jack-knife based on the average 
number of line-like features for all 150 GHz bolometers, 
as well as a jack-knife based on the bandwidth affected 
by the lines. Finally, we do an additional split using the 
average number of lines in an observation that appear to 
be related to the pulse-tube cooler. 

We calculate the \ 2 with respect to zero for each 
jack-knife over the range £ £ [2000, 10000] in bins with 
A£ — 500. Some of the tests are highly correlated. For 
example, we changed scan velocities approximately mid- 
way through the observations so splitting the data based 
on scan velocity is nearly identical to splitting the data 
between the first and second halves of the season. We 
calculate a correlation coefficient between the different 
tests by adding l/N b s for each common observation in 
a half, and subtracting l/N b s for each distinct obser- 
vation in a half. This algorithm returns unity for two 
identical sets and zero for two random sets, as we ex- 
pect 50% of the observations to be in common for two 
random selections. The correlation coefficients between 
the 13 jack-knives range from to 0.83, with the maxi- 
mum correlation being for the previously mentioned scan 
velocity and first half - second half jack-knives. We in- 
vert the jack-knife correlation matrix, C, and calculate 
X 2 = Via(C~ l )ijVj a . Here w iQ is the ratio of the band- 
power over the uncertainty for the i th jack-knife and a th 
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£-bin. The probability to exceed the measured % 2 for 
the complete set of thirteen jack-knives is 77% for the 
150 GHz data, 32% for the 220 GHz data and 57% for 
the combined set with both frequencies. We thus find no 
evidence for systematic contaminants in the SPT data 
set. 

5. FOREGROUNDS 

The main foregrounds at frequencies near 150 GHz are 
expected to be galactic dust emission, radio sources, and 
dusty star forming galaxies (DSFGs). Note, however, 
that the SPT field is selected to target one of the cleanest 
regions on the sky for g alactic dust emission, and in the 
Finkbeiner et al. ( 1999 1 model, dust emission is primarily 
on large angular scales. The contribution for the selected 
field on arcminute scales is insignificant. The primary 
foregrounds of consideration for this analysis are radio 
sources and DSFGs. 

Tens of bright radio sources are detected in the SPT 
maps at > 5 a, and contribute substantial amounts of 
power at both 150 and 220 GHz. Information on the 
fluxes and spectral indices of these and other sources 
significantly detected in the SPT maps can be found 
in V09. Without masking, the measured point source 
power is C7 imaskcd = 2.1 x 10"VK 2 at 150 GHz and 
^unmasked = lg x at 2 20 GHz. These esti- 

mates are dominated by the brightest few sources and 
thus subject to very large sample variance. We mask 
all sources with 150 GHz fluxes above the 5 a detection 
threshold, S — 6.4 mJy. By masking these sources we 
reduce the radio source contribution to the SPT band- 
powers by several or ders of magnitude. According to the 
de Zotti et al. (2005 1 model source counts, after masking 
bright sources, we expect a residual radio source contri- 
bution of C £ radio = 8.5 x 10" VK 2 at 150 GHz. 

The point source masking will remove the SZ contribu- 
tion from only a few galaxy clusters, leading to negligi- 
ble reduction of the SZ power. This is because the large 
majority of radio sources and DSFGs reside outside of 
SZ clusters. Most of the masked sources are identified 
as radio sources in V09. Extrapolations from lower fre- 
quency observations imply that, at 150 GHz, less than 
3% of clusters contain radi o source flux exceeding 20% 
of the tSZ flux decrement ( Lin et aL]|2009 Sehgal et al. 
20091. The masked sources are selected as increments 
at 150 GHz and therefore have fluxes much greater than 
20% of the tSZ of any associated galaxy cluster. The 
number of clusters masked by the radio source mask- 
ing should then be much less than 3% and negli gible. 
We also compar e the tSZ power spectrum in the Seh- 
gal et al. ( |2009[ ) simulated sky maps with and without 
masking > 6.4 mJy sources and find the difference to be 
<< 1%. A small number (six) of the masked sources 
are identified as DSFGs (V09). Galaxy clusters have a 
DSFG abundanc e only twenty times larger than the field 



(Bai et al. 20071, although they exceed the mass density 
of the field by a factor of 200 or more. Given the relative 
rarity of galaxy clusters, it follows that only a small frac- 
tion of DSFGs can live in galaxy clusters. Therefore, the 
number of clusters masked along with the DSFGs should 
be much smaller than six and negligible. 

Both the DSFG and radio source arguments above de- 
pend implicitly on the impact of p otentially masking 
^ 10 clusters. As a worst-case study, |Shaw et al.| ( |2009| 



consider the impact of masking the most massive ten 
clusters in the field and show that even in this extreme 
case, the tSZ power spectrum at £ = 3000 is reduced by 
only 11%. Of course, the point source masking will not 
select the most massive clusters and is highly unlikely to 
remove as many as ten clusters. Hence, the true impact 
will be significantly less. 

After we mask sources above the 5 a threshold, DSFGs 
are the dominant point source population in the SPT 
maps. These sources have been extensively studied a t 
higher fre quencies by SCUBA (|Holland et al.||1999| , 
MAM BO ( Kreysa et aI1|1998|, Boloca m ( [Glenn et al 
[19981), LA BOCA ( |Si7niIo^ral.|[2009t, A zTEC ( |ScW 
|et al. |2008|), SCU BA-2 polland et al.|2006l , and BEAST 
(JFascale et al. 20081, and there have been some pre 



liminary indications of their contributi on in previous 



small-scale power spectra at 150 GHz (Reichardt et al. 

|2009a|b[ ). 
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with the exact frequency dependence a function ot the 
dust emissivity, the dust temperature, and the redshift 
distribution of the galaxy population. This range of spec- 
tral indices corresponds to point source amplitude ratios, 
ST^q/ST^q, of 2.1-2.6 when expressed in units of CMB 
temperature. The measured spectral index of the DSFGs 
in the SPT maps is discussed extensively in H09. 

In order to obtain an unbiased estimate of the SZ power 
spectrum, it is essential that we take these sources into 
account in our fits and modeling. After masking the 
bright point sources, we significantly detect a Poisson dis- 
tributed power at 150 GHz of Cf = 7.1 ± 0.5 x 10~ 6 /iK 2 
(H09). This unclustered point source power climbs with 
increasing I to be comparable to the SZ effect by £ — 
2500 — 3000, and is the dominant astronomical signal in 
the maps at arcminute scales. 

The distribution of DSFGs on the sky is also expected 
to be clustered, resulting in a significant increase in 
power at £ < 3000. BLAST recently detected this clus- 
tered term for DS FGs at 600 - 1200 GHz (500 - 250 /xm) 
(Viero et al. |2009 1. Extrapolating the measured cluster- 
ing to 150 GHz, we expect the clustered contribution to 
be comparable to the tSZ effect. H09 have analyzed the 
SPT 150, 220, and 150 x 220 GHz bandpowers presented 
in $6l and have also found that the amplitude of the clus- 
tered component is indeed comparable to the tSZ am- 
plitude. Discriminating between clustered DSFGs and 
the tSZ effect would be extremely difficult for a single- 
frequency instrument as the angular dependencies are 
very similar. However, the two frequencies used in this 
analysis allow the spectral separation of the these two 
astronomical signals. 

6. RESULTS 

In addition to presenting the SPT angular power spec- 
tra, we wish to determine the amplitude of the tSZ 
power spectrum and to use this amplitude to constrain 
the normalization of the matter power spectrum, erg. 
This measurement requires separating the tSZ signal 
from the other astrophysical signals in our data which 
include primary CMB anisotropy (including lensing ef- 
fects), DSFGs (both Poisson and clustered components), 
and anisotropy due to the kSZ effect. The primary CMB 
anisotropy and Poisson point source component can be 
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Fig. 2. — Jack-knives for the SPT data set at 150 GHz (blue circles) and 220 GHz (black diamonds). For clarity, the 220 GHz 
jack-knives have been shifted to the right by A£ = 100. Top panel: Bandpowers of the "first half - second half jack-knife compared to the 
expected error bars about zero signal. Disagreement with zero would indicate cither a noise misestimate or a time-dependent systematic 
signal. Second panel: Power spectrum of the left-going minus right-going difference map. This test yields strong constraints on the accuracy 
of the detector transfer function deconvolution and on possible directional systematics. Third panel: Bandpowers for the difference map 
when the data is split based on azimuth. Signals fixed in azimuth such as side-lobe pickup from the nearby support building would produce 
non-zero power. We see no evidence for ground-based pickup. The cumulative probability to exceed the \ 2 observed in these three tests is 
76% at 150 GHz and 22% at 220 GHz. Bottom panel: The un-differenced SPT power spectra at each frequency for comparison. 



separated from a tSZ-like component on account of the 
distinct angular power spectra of these three signals. 
However the tSZ, kSZ, and clustered DSFG components 
are all expected to be roughly flat in T>£, and we de- 
pend on their distinct frequency dependences to separate 
them. We use a combination of the two frequency bands 
to remove the DSFG contribution to the power spectra. 
We address the remaining degeneracy between the tSZ 
and kSZ effects by repeating the analysis for a range of 
assu med kSZ models. 

In §6.1| we present bandpowers for 150, 2 20 a nd 150 x 
220 GHz SPT angular power spectra. In §6.2| we com- 
bine the two frequency bands to g enerate a set of DSFG- 
subtracted bandpowers. In |6.3| we describe the Monte 



Carlo Markov chain (MCMC) analysis used to estimate 
the tSZ power spectrum amplitude, parametrized as the 
normalization A^z of a mode l tem plate, from the DSFG- 
subtracted bandpowers. In ^ 6.4 we discuss the implica- 



tions of the measured tSZ power spectrum amplitude for 
ag and modeling of the tSZ effect. 

6.1. Power Spectra 

Figure [3] shows the bandpowers we compute by apply- 
ing the analysis methods described in |3|to one 100 deg 2 
field observed by SPT at 150 and 220 CIHz. The band- 
powers for the two frequencies and their cross-spectrum 
are tabulated in Table [TJ We combine these multi- 
frequency data as described in [j3]to produce the 'DSFG- 
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subtracted' bandpowers listed in Table |2] This power 
spectrum is compared to the results from WMAP5, 
ACBAR, and QUaD in Figure g] The best-fit model to 
this combined data set including the primary CMB, kSZ, 
tSZ, and a Poisson point source contribution is shown for 
reference. The primary CMB anisotropy is estimated for 
a spatially-flat, ACDM model, which includes gravita- 
tional lcnsing. The bandpower uncertainties are derived 
from the combination of simulations and the measured 
intrinsic variations within the SPT data as described in 
f|3] The bandpowers can be compa red to theo ry us- 
ing the associated window functions (Knox 1999). The 
bandpowers, uncertainties, and window functions may be 
downloaded from the SPT website. 25 

The SPT bandpowers ae dominated by the damping 
tail of the primary CMB anisotropy on angular multi- 
poles 2000 < I < 3000. At these multipoles, the band- 
powers are in excellent agreement with the predictions 
of a ACDM model determined from CMB observations 
on larger angular scales. On smaller scales, the SPT 
bandpowers provide new information on secondary CMB 
anisotropics and foregrounds which dominate the pri- 
mary CMB anisotropy. The SPT data presented here 
represent the first highly significant detection of power 
at these frequencies and angular scales where the pri- 
mary CMB anisotropy is sub-dominant. After masking 
bright point sources, the total signal-to-noise ratios on 
power in excess of the primary CMB are 55, 55, and 45 
at 150, 150 x 220, and 220 GHz respectively. The ma- 
jority of the high-€ power can be attributed to a Pois- 
son distribution of point sources (likely DSFGs) on the 
sky. H09 use the bandpowers presented in this work to 
measure both Poisson and clustered contributions to the 
measured powers and investigate the implications for the 
properties of DSFGs. The largest source of secondary 
CMB anisotropy at 150 GHz is expected to be the SZ 
effect, and we investigate SZ constraints in the following 
sections. 



6.2. DSFG- subtracted Bandpowers 

Our immediate goal is to measure the amplitude of the 
tSZ power spectrum. However, several signals in these 
maps have similar angular power spectrum shapes, and 
the single-frequency maps only constrain the sum of the 
power from these sources. For instance, the 150 GHz 
data effectively constrain the sum of the tSZ, kSZ, and 
clustered DSFG power. 

As discussed earlier, each of these components has a 
distinct frequency dependence so a linear combination of 
SPT's two frequency bands can be constructed (following 
{ 3.2.4 ) to minimize any one of them. H09 find significant 
evidence for a clustered DSFG power contribution to the 
single-frequency bandpowers listed in Table [T] with an 
amplitude comparable to that of the tSZ effect. We ex- 
pect the kSZ effect to be smaller than the tSZ effect at 
150 GHz on theoretical grounds. Additionally, due to the 
frequency dependencies of the components, removing the 
kSZ effect would inflate the relative contribution of clus- 
tered DSFGs with respect to the tSZ effect. Therefore, 
we choose to remove DSFGs from the SPT bandpowers. 

For a mean DSFG spectral index, a, the proper weight- 
ing ratio, x, to apply to the 220 GHz spectrum for DSFG 

25 http:/ /pole. uchicago.edu/public/data/lueker09/ 
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The spectrum in Table [2] and Figures [4] and [5] is pro- 
duced with a weighting factor of x — 0.325 correspond- 
ing to a mean spectral index of a = 3.6. The contri- 
bution from DSFGs can be completely removed only if 
every galaxy has the same spectral index. However, the 
comparative closeness of SPT's two frequency bands en- 
sures that power leakage into the difference maps remains 
small despite some expect ed sca tter in the spectral index 
of the dusty galaxies. In { 6.2.1 we motivate this choice 
of x and discuss predictions for the residual point source 
power. 

We assume that there is no tSZ contribution to the 
220 GHz data as the 220 GHz band is designed to be 
centered on the SZ null. Fourier transform spectroscopy 
measurements of the 220 GHz band pass confirm that 
the tSZ amplitude in the 220 GHz band will be < 5% of 
the 150 GHz amplitude. Any error incurred by subtract- 
ing roughly one third of the 220 GHz amplitude from 
the 150 GHz data would be less than 3%, far below the 
present ~ 40% statistical uncertainty on Agz (see Table 
3). 

It is important to note that the apparent tSZ power 
in the DSFG-subtracted bandpowers will be a factor of 
1/(1 - x) 2 = 2.2 higher than at 150 GHz as the differ- 
enced bandpowers have been normalized to preserve the 
amplitude of the primary CMB anisotropy. In this work, 
we report SZ amplitudes scaled to 153 GHz which is the 
effective band center of the 150 GHz band for a tSZ spec- 
trum. 

6.2.1. Residual Point Source Power 

The DSFG-subtracted maps have substantially less 
power due to both unclustered and clustered point 
sources as seen in Figure [5j However, we include a Pois- 
son point source amplitude in all fits since we expect a 
small fraction of the point source power to remain in the 
DSFG-subtracted maps. The best-fit amplitude is con- 
sistent with zero and unphysical negative values of the 
Poisson point source power are allowed due to noise. To 
prevent this, we place a prior on the residual Poisson 
point source power based on what we know about the 
observed DSFG populati on from H09 and rad io source 
population from V09 and de Zotti et al. (20051. 

The residual power in the subtracted map due to 
Poisson-distributed DSFGs, Cf SFG , depends on the scat- 
ter in spectral indices o~ a , the accuracy to which the 
mean spectral index a is known, and an estimate of the 
Poisson DSFG power in the 150 GHz band, Cf 5 ' 150 . For 
a given combination of these parameters, this residual 
DSFG power will be: 



^DSFG _ ^ps.150 



(26) 



°i [^(^150/^220)] + 1 



where v\^q and ^220 are the effective bandcenters of the 
150 and 220 GHz bands, and x and xtrac are the assumed 
and true values of map weighting ratio. 
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Fig. 3.— The SPT 150 GHz (purple circles), 150 X 220 GHz (green diamonds) and 220 GHz (blue triangles) bandpowers. The 
damping tail of the primary CMB anisotropy is apparent below t = 3000. Above I = 3000, there is a clear excess with an angular 
dependence consistent with point sources. These sources have low flux (as sources with > 6.4 mjy at 150 GHz have been masked) and 
a rising frequency spectrum, consistent with our expectations for Poisson distributed DSFGs. The point source population and resulting 
contributions to anisotropy power are discussed in more detail in H09. 



TABLE 1 
Single-frequency bandpowers 









150 GB 


z 


150 x 220 GHz 


220 GHz 


i range 


'off 


q(nK 2 ) a 


(pK 2 ) 


q ( M K 2 ) 


a ( M K 2 ) 


q ( M K 2 ) 


a (^K 2 ) 


2001 


- 2200 


2058 


240.5 


11.4 


269.5 


15.3 


336.2 


30.4 


2201 


- 2400 


2276 


139.4 


7.0 


155.0 


9.8 


195.4 


20.5 


2401 


- 2600 


2474 


114.3 


5.2 


119.2 


8.1 


197.6 


18.8 


2601 


- 2800 


2677 


79.0 


4.3 


114.0 


6.7 


213.3 


19.2 


2801 


- 3000 


2893 


54.1 


3.7 


77.7 


6.1 


183.1 


19.5 


3001 


- 3400 


3184 


47.0 


2.4 


76.5 


4.2 


183.7 


13.4 


3401 


- 3800 


3581 


36.9 


2.4 


75.8 


4.3 


207.7 


15.6 


3801 


- 4200 


3992 


35.0 


2.8 


82.7 


5.0 


276.5 


18.8 


4201 


- 4600 


4401 


33.9 


3.2 


87.5 


5.9 


252.8 


21.0 


4601 


- 5000 


4789 


33.6 


4.2 


99.4 


7.3 


306.1 


24.5 


5001 


- 5900 


5448 


47.1 


3.5 


135.1 


6.3 


433.6 


21.7 


5901 


- 6800 


6359 


60.5 


5.5 


158.2 


9.4 


469.6 


33.0 


6801 


- 7700 


7255 


93.5 


8.8 


195.8 


14.5 


579.7 


46.4 


7701 


- 8600 


8161 


76.8 


13.5 


221.8 


23.6 


776.5 


73.9 


8601 


- 9500 


9059 


115.6 


20.7 


299.7 


33.7 


986.1 


105.9 



Note. — Band multipole range and weighted value l e g, bandpower qs, and uncertainty cjb for the 150 GHz auto-spectrum, 
cross-spectrum, and 220 GHz auto-spectrum of the SPT field. The quoted uncertainties include instrumental noise and the 
Gaussian sample variance of the primary CMB and the point source foregrounds. The sample variance of the SZ effect, beam 
uncertainty, a nd ca libration uncertainty is not included. Beam uncertainties are shown in Figure [I] and calibration uncertainties 
are quoted in §2.2| Point sources above 6.4 mjy at 150 GHz have been masked out in this analysis. This flux cut substantially 
reduces the contribution of radio sources to the bandpowers, although DSFGs below this threshold contribute significantly to 
the bandpowers. 



We examine the residual Poisson point source ampli- 
tudes for a broad range of weighting ratios to estimate 
the optimal x value and the error in that estimate. For 
each x, we estimate the probability that Cf s (x) is less 
than Cf s (x — 0.325) using the MCMC chains described 



in |6.3| The resulting probability distribution is taken 
to be the likelihood function for x trU c- As shown in Fig- 
ure |6l there is a broad maximum for x — 0.25 to 0.4 and 
we adopt the best fit, x = 0.325, for the following results. 
In the absence of a direct measurement, we place a con- 
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TABLE 2 

DSFG-SUBTRACTED BANDPOWERS 



t range 


7 

tefl 


q {[xK ) 




2001 - 


2200 


2058 


221.3 


16.9 


2201 - 


2400 


2276 


130.2 


11.2 


2401 - 


2600 


2474 


126.5 


10.3 


2601 - 


2800 


2677 


60.2 


7.7 


2801 - 


3000 


2893 


50.4 


8.0 


3001 - 


3400 


3184 


36.6 


5.9 


3401 - 


3800 


3581 


21.0 


6.5 


3801 - 


4200 


3992 


22.9 


8.4 


4201 - 


4600 


4401 


8.1 


9.5 


4601 - 


5000 


4789 


3.0 


11.3 


5001 - 


5900 


5448 


11.2 


10.5 


5901 - 


6800 


6359 


16.0 


16.2 


6801 - 


7700 


7255 


60.3 


27.9 


7701 - 


8600 


8161 


32.1 


42.9 


8601 - 


9500 


9059 


54.7 


63.8 



Note. — Band multipole range and weighted value £ c ft, bandpower qs, a nd unc ertainty ob for the DSFG-subtracted maps 



of the SPT field. These bandpowers correspond to a linear combination (see [ 3.2.4 1 of the 150, 150 x 220, and 220 GHz power 
spectra, optimized to remove emission from DSFGs below the point source detection threshold of SPT. Point sources above 
6.4 mjy at 150 GHz have been masked out in this analysis. The quoted uncertainties include instrumental noise and Gaussian 
sample variance of the primary CMB and point source foregrounds. The sample variance o f the SZ e ffect , beam uncertainty and 
calibration uncertainty is not included. Beam and calibration uncertainties are quoted in §2.1| and §2.2| and shown in Figure [5] 



servative uniform prior on the scatter in DSFG spectral 
indices, 0.2 < a a < 0.7, as discussed in H09. 

Our expectation for the residual radio contribution to 
the DSFG-subtracted b andpowers, C^ adl °, is based on 
the de Zotti et al.| ( p005[ | radio source count model. This 
model is in excellent agreement on the high flux end with 
the SPT source counts (V09). As discussed in H09, the 
residual radio source power after masking is expected to 
be a small fraction of the DSFG power at 150 GHz. How- 
ever, this small radio contribution may be comparable to 
the residual DSFG power in the DSFG-subtracted spec- 
trum. The radio source power c an be calculated from 
the integral of S 2 dN/dS for the jde Zotti et al.| ( |2005 1 
counts model from zero to the flux masking threshold 
of 6.4 mjy. We compute the power these sources con- 
tribute to the optimal DSFG-subtracted spectrum to be 
3.9 x 10 -7 /xK 2 by assuming an average spectral index of 
a = —0.5 based on the detected sources in V09. This 
power level is near ly identical to that predicted by the 
Sehgal et al. (20091 simulations. To allow for a variation 
in spectral index as well as uncertainty in the model nor- 
malization when extended to lower flux sources, we as- 
sign a conservative uncertainty of 50% on the predicted 
residual radio source power in the prior. 

We combine this information to create a prior on 
the residual point source power in the DSFG-subtracted 
maps Cg S = Cf + C^ adl °. This prior spans the range 
Cf e[3.5, 9.0]xlO-VK 2 at 68% confidence and [1.2, 
13.9] xlO-VK 2 at 95% confidence. The best-fit value 
of the residual Poisson component before applying the 
prior is C^no prior) = (6.2±6.4) x 10~ 7 /iK 2 and lies at 
the middle of our assumed prior range. The upper end 
of the 95% range is approximately 20% of the best-fit 
value of the Poisson point source power in the undiffcr- 
enced 150 GHz bandpowers. This suggests that we have 
subtracted over 80% of the point source power from the 



150 GHz spectrum, with the residual point source power 
largely from radio sources. Without this prior on the 
Poisson point source amplitude, the uncertainty on the 
Asz detection presented in the next section would in- 
crease by ~ 50%. 

6.2.2. Residual Clustered Point Source Power 

We assume that the contribution of clustered point 
sources is insignificant in the DSFG-subtracted band- 
powers. Using a combination of the SPT bandpowers 
at 150, 220, and 150 x 220 GHz as well as the DSFG- 
subtracted spectrum, H09 argue that the residual clus- 
tered DSFG component in the DSFG-subtracted band- 
powers is less than 0.3 /iK 2 at 95% confidence. This is 
several percent of the SZ power spectrum but negligible 
at the current detection significance of < 3 a. We also 
argue in Sj5] that clustered radio sources are negligible. 
Therefore, residual power from clustered point sources 
will not bias SZ constraints from the DSFG-subtracted 
bandpowers. 

The above argument holds if the point sources are un- 
correlated with the SZ signal. However, if the clustered 
term was completely anti-correlated with the SZ signal, 
the measured SZ power in Table [3] could underestimate 
the true SZ power by 38%. This is unlikely for two rea- 
sons. First, the residual after DSFG subtraction should 
be uncorrelated as long as the spectral dependence of 
cluster member DSFGs is similar to the general DSFG 
population. Second, as argued in f|5] most DSFGs are 
not galaxy cluster members. We also lo ok at the corre- 
lations between the |Sehgal et al~| (|2009|) DSFG a nd tSZ 
simulated sky maps The |Sehgal et ai.| ( |2009| ) DSFG 
model scales the number density ot DSFG cluster mem- 
bers linearly with cluster mass ; this is a subs tantially 
stronger scaling than observed (|Bai et al.||2007 1. There- 
fore, estimating the cluster-DSFG correlation from the 
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Fig. 4. — WMAP5 (blue squares), ACBAR (green triangles), QUaD (turquoise diamonds) and the SPT (black circles) DSFG- 
subtracted SPT bandpowers are plotted over the best-fit, lensed ACDM cosmological model (dashed red line), best-fit tSZ power 
spectrum (solid black line), homogeneous kSZ model (dashed black line), and residual Poisson-distributed point source contribution 
(solid orange line). The combined best-fit model is shown by the solid red line. The plotted SPT bandpowers have been multiplied 
by the best-fit calibration factor of 0.92. Point sources above 6.4 mjy at 150 GHz have been masked. The patchy kSZ template is also 
shown for reference (dotted black line). The DSFG-subtracted bandpowers are normalized to preserve the amplitude of the primary 
CMB anisotropies. 



Sehgal et al. (20091 should be overly conservative. We 
calculate the anti-correlation coefficient between the ISe-l 



hgal et al. (2009) simulated tSZ maps at 148 GHz and 
a linear combination of the 148 and 219 GHz simulated 
IR source maps with the same weighting as used for the 
SPT DSFG-subtracted bandpowers. We find the anti- 
correlation coefficient between the tSZ effect and total 
DSFG power to be 21% using the power measured at 
£ = 3000. This highly conservative upper limit of 21% 
on the anti-correlation implies that the true SZ power is 
underestimated by less than 4%. Based on these argu- 
ments, we assume that correlations between SZ signals 
and emission from cluster member galaxies is negligible 
in this analysis. 

Separately, one might worry about correlations be- 
tween radio sources and SZ clusters. Radio sources could 
suppress the tSZ signal by "filling in" the t SZ d ecre- 
ments. However , the work by Lin et al. (2009) and |Seh- 
gal et al. (20091 shows the number of cluster-correlated 
radio sources is expected to be small. For instance, we 
can examine this correlation in the simulate d tSZ and 

^009" 



radio source maps produced by Sehgal et al 
look at the power spectrum at £ 



We 

MKTfcFlhe tSZ 



map, the tSZ+radio-source map, and the radio source 
only map after masking sources above 6.4 mjy. We find 
an anti-correlation coefficient of 2.3% for the two com- 
ponents. Given the expecte d rad io source power level 
of ~ 0.6 ^K 2 at t = 3000 (j |6 2.l| ), radio source-cluster 
correlations should not affect the results in this work. 

6.3. Markov Chain Analysis 



The DSFG-subtracted bandpowers presented in {6.2 
detect at high significance a combination of the primary 
CMB anisotropy, secondary SZ anisotropies, and residual 
point sources. In this section, we use an MCMC analysis 
to separate these three components and to produce an 
unbiased measurement of the tSZ power spectrum am- 
plitude. 

6.3.1. Elements of the MCMC Analysis 

We fit the DSFG-subtracted bandpowers to a model in- 
cluding the lensed primary CMB anisotropy, secondary 
tSZ and kSZ anisotropies, and a residual Poisson point 
source term. We use the standard, six-parameter, spa- 
tially flat, lensed ACDM cosmological model to predict 
the primary CMB temperature anisotropy. The six pa- 
rameters are the baryon density Qf,, the density of cold 
dark matter f2 c , the optical depth to recombination t, 
the angular scale of the peaks 0, the amplitude of the 
primordial density fluctuations ln[10 1( M s ], and the scalar 
spectral index n s . To fit the high-i? power, we extend the 
basic, six-parameter model with two additional parame- 
ters: the amplitude of a tSZ power spectrum template, 
Asz, and a constant, Cf s , such as would be produced by 
a Poisson distribution of point sources on the sky. We 
also explore the potential impact of the kSZ effect on 
these parameters by using three different kSZ models. 

Gravitational lensing of CMB anisotropy by large scale 
structure tends to increase the power at small angular 
scales, with the potential to influence a SZ power spec- 
trum measurement. The calculation of lensed CMB spec- 
tra out to I = 10000 proved prohibitively expensive in 
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Fig. 5. — The SPT 150 GHz (purple diamonds) and DSFG-subtracted (black circles) bandpowers over-plotted on the best-fit models 
to the DSFG-subtracted bandpowers. The best-fit, lensed ACDM cosmological model for the primary CMB anisotropy is shown by the 
dashed red line, while the sum of the best-fit ACDM model, kSZ, tSZ and point source terms is shown by the solid red line. The 
primary CMB anisotropy alone is a poor fit to the SPT data. The uncertainties on the DSFG-subtracted bandpowers are larger for two 
reasons. First, the normalization convention inflates the uncertainties by a factor of 1/0.675 2 , and second, these bandpowers also include the 
more noisy 220 GHz data. Beam and calibration uncertainties are marked by a second blue error bar for the DSFG-subtracted bandpowers 
only. Note that the calibration and beam uncertainties are correlated between ^-bins. The 150 GHz data has been shifted to the right by 
A£ = 40 for clarity. Point sources above 6.4 mjy at 150 GHz have been masked in this analysis. 

spectrum for the same best-fit WMAP5 cosmology aver- 
ages ~ 8.6 /jK 2 near £ = 3000 where SPT has the high- 
est S/N on the SZ spectrum. In the diffe renc ed spectra 
used to derive the A$z constraints (see [6.2), the ratio 
of lensing to the tSZ effect is suppressed by a factor of 
(1 — x) 2 — 0.46. This reduction occurs because the sub- 
tracted spectrum is normalized such that the CMB power 
is unaffected, though the SZ power is enhanced. As a 
result, we expect possible lensing misestimates to intro- 
duce a negligible error of < 3% on the Asz constraints. 
Of course, this error will be larger for smaller values of 

The tSZ templat e we use as a fiducia l model is based 
on simulations by |Sehgal et al. (2009). 26 The simula- 
tions are for a WMAP5 cosmology with a% = 0.80 and 
fl^h = 0.0312 and an observing frequency of 148 GHz. 
We rescale the template to 153 GHz, the effective center 
frequency of the SPT 150 GHz band for signals with a 
tSZ spectrum. 

For a limited number of cases, we also compare the re- 
sults of the Sehgal numerical tSZ template to those of the 
Komatsu & Seljak analytic template with WMAP5 pa- 
rameters 27 j Komatsu fc Sel.iak ||2002|) and the numerical 




Fig. 6. — Probability that the residual point source power in the 
DSFG-subtracted map constructed by rhiso ~ xfft220 ls lower than 
the value at x = 0.325 as a function of x. We can interpret this as 
the probability that a given value of x is the true value. There is a 
broad maximum centered at 0.325, which corresponds to a spectral 
index of 3.6 between the 150 and 220 GHz bands. This is consistent 
with the ratio of point source power between single-frequency fits 
to the 150, 150 X 220, and 220 GHz bandpowers. We estimate the 
x which minimizes residual point source power to be 0.325 d_ 0.08. 

computational time. We avoid this computational limi- 
tation by calculating the lensing contribution for the the 
best-fit cosmological model, and adding this estimated 
lensing contribution to the unlensed CMB power spec- 
trum calculated at each step in the Markov chain. Given 
the small allowed range of f_ m with current CMB data, 
we predict that using the fixed lensing contribution will 
misestimate the actual lensing by less than 30%. We 
have checked this assertion on a sampling of parameter 
sets drawn from the chains. The lensing contribution to 
the high-€ spectrum is ~1.5 fiK 2 while the modeled tSZ 



template by phaw et al.| ( |2009[ ). The primary difference 
between the Shaw and bengal simulat ions is the value of 
the energy feedback parameter in the Bode et al. ( 2007 1 
intracluster gas model. The Shaw simulations have the 
feedback parameter reduced to 60% of the Sehgal value. 
Increasing the feedback parameter causes the gas distri- 



26 http:/ /lambda. gsfc.nasa.gov/toolbox/tb_cmbsim_ov.cfm 

27 http://lambda.gsfc.nasa.gov/product/map/dr3/pow_sz_spec_get.cfm 
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bution in clusters to 'puff-out' or inflate; in low mass 
clusters the gas may become unbound altogether. The 
overall effect is to reduce the predicted tSZ power spec- 
trum, especially at small angular scales. We limit our 
analysis of specific tSZ models to these three models. 

An independent measurement of the kSZ spectrum is 
outside the scope of this work, however, it is necessary to 
take the kSZ effect into account when inferring the tSZ 
amplitude from the data. We consider three kSZ cases 
based on two published kSZ models. The kSZ effect is 
assumed to be zero in the first case ("no kSZ" ). As a sec- 
ond, in termediate case, we use the model by |Sehgal et al.| 
( 2009 1 (referred to as "homogeneous kSZ" ) . This model 
includes kSZ contributions from a homogeneous reion- 
ization scenario, but does not include the additional kSZ 
power produced by patchy reionization scenarios. We 
take this case to be the fiducial kSZ model. We include 
an estimate of patchy reionization in the third kSZ case. 
For the patchy reio nization phase, we use the "brief his- 
tory" model B from Zahn et al. (20051, re-calculated for 
WMAP-5 best fit cosmological parameters. The sum of 
the homogeneous kSZ model and the patchy contribution 
will be referred to as the "patchy kSZ" model. We do 
not scale these templates for different cosmological pa- 
rameter sets as we expect the kSZ theoretical uncertainty 
to be at least as large as the cosmological dependence. 
Finally, we include a residual Poisson point source com- 
ponent in all chains. The c onstruction of the point source 
prior is outlined in §6.2.1| 

Previous CMB experiments have produced exceptional 
constraints on the primary CMB anisotropy, and we us e 
the band powers from WMAP5 (|Dmrkley et al.||2009|), 

Brown 



2009a I, and QUaD 



ACBAR (|Reichardt et al., 

et al.|2009| at t < 2200 m all parameter fitting. We refer 
to this collection along with the SPT DSFG-subtracted 
bandpowers as the 'CMBalF data set. It is important 
to note that the two-parameter extension to the ACDM 
model for point sources and the tSZ effect is only re- 
stricted to the SPT bandpowers. This restriction is im- 
posed for two reasons. First, the point source contri- 
butions to each experiment may be different due to the 
different frequencies and flux cuts for masking sources. 
Arguably, the point source power is likely to be similar in 
the 150 GHz SPT, ACBAR, and QUaD results, but we 
would be unable to use frequency information to discrim- 
inate between the SZ effect and clustered DSFGs. Sec- 
ond, the primary CMB is dominant below I ~ 3000 and 
the other experiments lack sufficient statistical weight at 
high-^ to improve upon SPT's measurement of SZ effect 
and point source power. 

Parameter estimation is performed by MCMC sam- 
pling of the full multi-dimensional parameter space us- 



ing an exte nsion of the CosmoMC package (Lewis & 
Bridle 2002). We include t he code extension p roduced 



by the QUaD collaboration (Brown et al. 20091 to han- 
dle uncertainties on a non-Gaussian beam in CosmoMC. 
CMB power spectra for a given parameter set are cal- 
culated with CAMB ( iLewis et ah] |2000| ). We use 
the WMAP5 likelihood code publicly available from 
|http: / /lambda. gsfc .nasa.gov| After the burn-in period, 
each set of four chains is run until the largest eigen- 
value of the Gelman-Rubin test is smaller than 0.0005. 
Wide uniform priors are used on all six parameters of the 
ACDM model. A weak prior on the age of the Universe 



(to e [10,20] Gyrs) and Hubble constant (h G [0.4,1]) is 
included in all chains, but should not affect the results. 
We use a uniform prior on A§z over a wide range from 
-1 to 10 times the value expected for erg = 0.80. 

6.3.2. Constraints on SZ amplitude 

We fit for the normalization factor of a fixed tSZ tem- 
plate and, in Table [3j report both this template-specific 
normalization, A§z, and the total inferred SZ-power at 
£ = 3000, near the multipole with maximum tSZ detec- 
tion significance. This estimate of the SZ power includes 
both tSZ and kSZ terms, and is included to facilitate 
comparison with other SZ models. We expect both the 
thermal and the kinetic SZ spectra to vary slowly with 
angular multipole. 

The chains are run for three different assumptions 
about the kSZ effect: the no kSZ, homogeneous kSZ, 
and patchy kSZ models described above. In each case, 
we produce MCMC chains with a fixed kSZ amplitude. 
The x 2 °f the 15 SPT bandpowers is between 16.3 and 
16.4 for the best fits of the three kSZ cases considered; 
there is essentially no impact on the quality of the fit. 
The Ag-z, and power constraints for each case are listed 
in Table [3] and plotted in Figure [7] 

The bandpower uncertainties in Table [2]do not include 
the sample variance of the tSZ effect, and we must con- 
volve the Agz distribution in the chains with an esti- 
mate of the sample variance to find the true A$z like- 
lihood function. The sample varianc e of the tSZ effect 
is estimated using the simulations of Shaw et al. (2009) 
re run with the same in tracluster gas model parameters 
as Sehgal et al. ( 2009 1 . The simulation consists of 300 
map realizations of the size of the SPT sky patch. These 
are constructed from a base sample of 40 independent 
maps by separating the components of each map into 
eight redshift bins (between < z < 3) and shuffling 
these bins between maps to generate a larger sample. We 
take the spectrum of each realization and find the best- 
fit Agz amplitude after weighting £-b'ms based on the 
SPT bandpower uncertainties. We use the distribution 
of these amplitudes to map out the likelihood function 
for the tSZ power (see Figure [8| . Due to the number 
of independent realizations, wenave limited ability to 
resolve the tail of the likelihood function. Fortunately, 
the non-Gaussianity is small for such a large sky area 
and the distribution of ln(^lsz) is well-fit by a Gaussian 
with a 12% width. We use this Gaussian fit as an esti- 
mate of the full likelihood surface. Small deviations from 
the true description of the tSZ sample variance will not 
impact the final results as the sample variance is small 
compared to both the statistical uncertainties and the 
assumed 50% model uncertainty. The uncertainties on 
Asz are essentially unchanged by the inclusion of the 
tSZ sample variance. The sample variance and model 
uncertainty are shown in Figure [8j 

As discussed earlier, the DSFG-subtracted bandpowers 
are sensitive to a linear combination of the tSZ and kSZ 
effects. We expect analysis of 2009 and later SPT data 
which include 95 GHz data to be able to separate the 
two SZ effects. The linear combination is not a simple 
sum, as the frequency-differencing used to produce the 
DSFG-subtracted spectrum suppresses the kSZ relative 
to the tSZ by a factor of (1 - x) 2 = 0.46. This factor 
is uncertain at the 15% level due to the relative cali- 
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TABLE 3 



Constraints on Ag Z and <j$ 




primary CMB 


CMB ± A sz : 


A sz : 

Asz (w homogeneous kSZ): 
A sz (w patchy kSZ): 


- 


0.55 ±0.21 
0.42 ±0.21 
0.34 ±0.21 


SZ power at = 3000 
( tSZ + 0.46 x kSZ): 




4.2 ± 1.5 nK 2 


kSZ power at £ = 3000 
homogeneous kSZ 
patchy kSZ 




2.0 ^K 2 
3.3 /xK 2 


cr 8 (no kSZ): 

erg (w homogeneous kSZ): 
erg (w patchy kSZ): 


0.795 ± 0.033 
0.794 ±0.028 
0.788 ±0.029 


0.778 ± 0.024 
0.773 ±0.025 
0.770 ±0.024 



Note. — The 1 a constraints on <rg derived from the DSFG-subtracted analysis of the SPT data, when using the simulations 

in 



Shaw et al. (12009) to estimate the non-Gaussian cosmic variance of the tSZ power spectrum. The best-fit value for the 
amplitude of the tbZ power spectrum is also shown, normalized to unity for a WMAP5 cosmology with og = 0.8. Asz = 1 
corresponds to a power of 7.5 fjK 2 &t£ = 3000. Results are shown with cosmic variance added in quadrature to the statistical 
uncertainty, however, the Asz constrai nts are dominated b y statistical uncertainties. Results are shown for no kSZ effect, for a 
homogeneous model of the kSZ e ffect (jSehgal et al.||2009j ) and the homogeneous model with an additional patchy reionization 
power contribution (Zah n et al. | |2005p T Finally, the joint constraint on the combined kSZ/tSZ power is shown under the 
assumption that the two templates are effectively degenerate. For reference, we also quote the power of the two kSZ models 
considered. The SPT data constrains the combined amplitude of the SZ contributions, which we quote at £ = 3000 where the 
measurements have the most constraining power. The kSZ and tSZ receive different pre-factors in the frequency-differenced 
analysis due to their relative spectral dependence. 



bration uncertainty between the bands. The SPT data 
detect the combined SZ effect at 2.6 a with tSZ ± 0.46 
x kSZ = 4.2 ± 1.5 fiK 2 at t = 3000. Using this com- 
bined constraint implicitly assumes that the tSZ and kSZ 
templates are perfectly degenerate, which is a good as- 
sumption for the current data quality and templates used 
in this work. 

We can compare the power detect ed with SPT to tha t 
reported by the CBI collaboration (Sievers et al. 2009). 
We use the best-fit normalization of a Komatsu &; Sel- 
jak template for WMAP5 parameters to compare directly 
the results of the two experiments. There are sub-percent 
differences in the assumed values of a% and f2f,/i between 
the template we adopted here and that used in |Sievers 
et al.| (20091, which would change the amplitude of the 
template by ~1% for an assumed scaling of al(Qt,h) 2 . 
This effect is negligible. We find the best-fit normal- 
ization of the WMAP5 Komatsu & Seljak model to be 
0.37 ± 0.17 for the SPT data under the homogeneous 
kSZ scenario. This is 2.4 a below the best-fit CBI nor- 
malization of 3.5±1.3. The model includes the frequency 
dependence of the tSZ effect. The smaller SPT bandpow- 
ers suggest that the CBI excess power may be produced 
by foregrounds with a frequency dependence falling more 
steeply than the SZ effect such as radio sources. 

6.4. Implications of the Asz Measurement 

The best-fit normalization for the fiducial tSZ spec- 
trum, Agz, is significantly lower than unity. The cosmo- 
logical parameters assumed when generating this tem- 
plate may be slightly different than the best-fit models, 
so we scale the template to match the cosmologies ex- 
plored by the Markov chain. When these scalings are 
taken into account, the measured values of Asz are still 



low. The tSZ power spectrum depends on the details of 
how the baryon intracluster gas populates dark matter 
halos and on cosmology through the number density of 
these halos. The paucity of tSZ power may reflect an 
overestimate of the intracluster gas pressure by the fidu- 
cial model and so we compare the template to other SZ 
models. At the same time, this low value of Asz favors a 
shift in the derived cosmological parameters, particularly 
cr 8 . 

Even a ~ 2.5 a detection of Asz will produce cosmolog- 
ically interesting constraints on as, as we expect Asz to 
scale strongly with er 8 and less strongly with the baryon 
density. The scaling is approximately aj(flbh) 2 where 



7 < 7 < 9, depending on the exact cosmolo gy (Komatsu 



fc Kitayama|1999 Komatsu & Seljak 



2002). We explore 



the erg dependence under the Press-Scnechter halo model 
and find that this relationship steepens with the cur- 
rently favored lower values of <t 8 . Sampling cosmolog- 
ical parameter values from the WMAP5 MCMC chains 
(which properly treats degeneracies between c 8 and other 
parameters), we find the modeled amplitude of the tSZ 
power spectrum varies approximately as a\ l . Constrain- 
ing the amplitude of the tSZ effect offers an independent 
measurement of as that can be compared to measure- 
ments based on primary CMB anisotropy or large scale 
structure. Such comparisons test our understanding of 
the physical processes involved in structure formation. 

For each point in the MCMC chain, we calculate the 
predicted Asz value from the six basic cosmological pa- 
rameters. For this calculation, we use the mass function 
of 
ga 



Jenkins et al. (20011 to determine the abundance of 
axy clusters of a give n mass. We then use the mass- 
concentration relation of Duffy et al. (2008) to determine 
the dark matter halo properties, and the gas model used 
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Fig. 7. — The ID marginalized Agz constraints from the SPT DSFG-subtracted bandpowers. Three kSZ cases are considered: no kSZ 
effect (dashed line), the homogeneous kSZ mode l (solid line) and the homogeneous model plus a patchy reionization term (dotted line). 
These models are described more fully in j )6.3.1| Top axis: The corresponding tSZ power at i = 3000 for reference. The no-kSZ curve 
(dashed line) can be interpreted as a constraint on the sum of £>| § + 0.46 x D!? sz . 



in Komatsu & Seljak (20021 to estimate the tSZ signal 
for each halo according to its mass. In order to con- 
vert this analytic tSZ spectrum into an amplitude, we 
take an ^-weighted average designed to match the rela- 
tive weights each multipole receives in the real tSZ fits. 
This amplitude is normalized to unity for the cosmolog- 
ical parameters assumed in the fiducial tSZ model. We 
also allowed the mass-concentration index to vary with 
cosmology by appropriate scaling of the characteristic 
mass M 4) but this was found to be a negligible effect 
within the explored range in parameters. 

At each point in the chain, the measured tSZ ampli- 
tude is compared to the predicted tSZ amplitude to con- 
struct a tSZ scaling factor, Asz/Ag2° ory . This procedure 
will account for any correlations between the measured 
Asz parameter and the six ACDM parameters in a self- 
consistent fashion, although we do not see evidence for 
such correlations in the current data. The distribution of 
scaling factors vs. eg is illustrated by the black contours 
in Figure [9] 

In general, the tSZ scaling factors are less than unity. 
These low scaling factors suggest either an over-estimate 
of the tSZ effect or lower values of erg. Models predicting 
larger kSZ or tSZ effects lead to lower scaling factors. 
However, the results can not be purely explained by an 
over-estimate of the kSZ effect since this tension persists 
in the no-kSZ case. Alternatively, the scaling factors may 
indicate that the explored range in cosmological param- 
eters is systematically overestimating the RMS of the 
mass distribution. For instance, as we see in Figure [9j 



points of the chain with lower values of og, have scaling 
factors closer to unity. 

The first interpretation of the low tSZ scaling factor 
is that the Sehgal tSZ template overestimates the tSZ 
power spectrum. There is currently some degree of un- 
certainty in the expected shape and amplitude of the 
tSZ power spectrum as predicted by analytic models or 
hydrodynamical simulations. One reason for this is that 
cosmological simulations of the intracluster medium have 
only recently begun to investigate in detail the impact of 
radiative cooling, non-gravitational heating sources (such 
as AGN), and possible regulatory mechanisms between 
them. The computational expense of running hydrody- 
namical simulations with sufficient resolution to resolve 
small-scale processes (such as star-formation) while en- 
compassing a large enough volume to adequately sample 
the halo mass function is prohibitive to a detailed anal- 
ysis of the predicted SZ power spectrum. In order to ac- 
curately predict the tSZ power spectrum, it is especially 
important to correctly model the gas temperature and 
density distribution in low mass (M < 2 x 10 14 h~ 1 M Q ) 
and high redshift (z > 1) clusters, which contribute sig- 
nificantly to the power spect rum at the angular scal es 
where SPT is most sensitive (Komatsu & Seljak 2002). 

In Figure flO| we plot the tSZ power spectrum derived 
from severalaifferent simulations (note that all curves 
have been normalized to the fiducial cosmology). The 
thick black solid line sh ows the base templa te, obtained 
from maps generated by Sehgal et al.| ( |2009[ ) . The black 
dot-dashed line shows the power spectrum obtained from 
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Fig. 8. — Sample variance and assumed theoretical uncertainty on the tSZ amplitude. The histogram shows the sample variance of 
^Sz/(^4sz)i where (Agz) is the mean value measured over a sample of 300 simulated maps. The overlying dashed black line shows the 
lognormal fit to the distribution, i.e., a Gaussian fit to ln(j4sz)), with <j\ u a = 0.12. The solid black line is the assumed 50% theoretical 
uncertainty on Agz, which we model as a Gaussian distribution in ln(Agz). 



the simulations of Shaw et al. (2009 1 , which have a lower 
energy feedback parameter than the Sehgal ct al.| ( |2009| ) 
simulations. Increasing the amount of feedback energy 
has the effect of inflating the gas distribution and sup- 
pressing power at small angular scales. The red solid and 
dot-dashed lines show the power spectrum from maps 
constructed from a 240/i _1 Mpc box simul ation run using 



mogeneous kSZ model, and 0.34 ± 0.21 for the patchy 
kSZ model. The tensi on grows worse fo r the medium- 
feedback simulatio ns bylShaw et al.l (12009 1. Although the 



the E ulcrian hydrodynamics code CART ( jKravtsov et al. 
|2005| Douglas Rudd, private communicat ion) a nd from 



the MareNostrum simulation ( Zahn et al.||2010"| respec- 
tively. Both of these simulations were run m the adia- 
batic regime, i.e., no cooling, star-formation or feedback. 
The gas distribution is more centrally concentrated than 
in the Shaw or Sehgal models, producing tSZ power spec- 
tra with significantly more power at small angular scales 
and less at larger angular scales. The blue solid line 
shows the analy tic template predicted by the |Komatsu| 
fc Seljak] ( |2002[) halo model calculation. The histogram 
shows the SPT sensitivity to the thermal SZ signal in 
each ^-band (with arbitrary normalization). 

The current data are in tension with even the high- 
feedback simulations for the CMB-derived best-fit cos- 
mological parameter set. Even with the kSZ effect set to 
zero, the tSZ scaling factor is only 0.55 ± 0.21 of what 
is predicted for the fiducial WMAP5 cosmology. Mean- 
while, the tSZ scaling factors are 0.42 ± 0.21 for the ho- 



Bode et al. (20071 model is calibrated to reproduce ob- 
served x-ray scaling relations for high-mass, low-redshift 
clusters, it may significantly over-estimate the contribu- 
tion of low-mass or high-redshift clusters (for which there 
are few direct x-ray observations with which to compare) . 

As mentioned above, an alternate interpretation of the 
low tSZ scaling factor is that the CMBall parameter 
chains without SZ constraints overestimate <r$. In or- 
der to constrain <7g based on the measurement of Agz, 
we need to construct the likelihood of observing Agz for 
a given cosmological parameter set. At a minimum, this 
likelihood would reflect the 12% uncertainty due to sam- 
ple variance of the tSZ effect. However given the large 
range of tSZ predictions, we add an additional theory un- 
certainty in quadratur e w ith the sample variance. The 
grey region in Figure [lO] shows the 1 a region encom- 
passed by a lognormal distribution of width aAsz = 0.5 
around the fiducial model (black solid line) . In the range 
where SPT is most sensitive (2000 < t < 6000), all the 
predicted curves lie within this region. In order to ac- 
count for the range of predicted power spectra, we adopt 
a 50% theory uncertainty on the value of ln(A sz ) in the 
likelihood calculation. The resulting a% constraints are 
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Fig. 9. — Two-dimensional likelihood contours at 68% and 95% confidence for ag versus the tSZ scaling factor, Agz/ A l c^ ory , derived 
from the SPT DSFG-subtracted bandpowers. For each point in the Markov chain, the tSZ scaling factor compares the Agz value fit to the 
SPT data to the A l ^° Ty vame predicted for that point's ACDM model parameters (see [6.4 1. The black contours show the likelihood 
surface for the CMBall dataset. We observe no dependence between Agz an d the six parameters of the ACDM model. The tilt towards 

higher scaling factors at lower ag is expected since the predicted ^sz GOry depends steeply on the value of <rg. The black contours also do 
not account for the cosmic variance of Agz't without cosmic variance or uncertainty in modeling the tSZ power, the tSZ scaling factor would 
be constrained to be exactly unity (solid orange line). The red shaded regions about unity illustrate the uncertainty we assume for 
the tSZ scaling factor due to theoretical uncertainty and cosmic variance (also see Fig. 111. This uncertainty is modeled as a log- normal 
distribution. The measured value of the tSZ scaling relation, including the theoretical uncertainty and sample variance in the model, is used 
to importance sample the Markov chain and obtain the likelihood surface marked by the blue contours. Left panel: Likelihood surfaces 
assuming no kSZ contribution. Right panel: Likelihood surfaces assuming the patchy kSZ model. The constraints for the homogeneous 
kSZ model will lie between the results for these two cases. 



dominated by the theory uncertainty. With this prior, we 
construct a new chain from the original parameter space 
through importance sampling. The regions preferred by 
this prior are shown in Figure [9] as are the results of the 
new Markov chain. 

Constraints on as with and without including the S PT 
Asz measurements are shown in Tableland Figure 11 



Under the assumption of the homogeneous kSZ model 
and a 50% theoretical uncertainty in the amplitude of the 
tSZ powerspectrum, the addition of the SPT data slightly 
tightens the constraint on as, while reducing the central 
value from a 8 = 0.794±0.028 to 0.773±0.025. The uncer- 
tainty in the resulting constraint on as is dominated by 
the large theoretical uncertainty in the tSZ amplitude. If 
instead we assume that the fiducial tSZ model is perfectly 
accurate and do not account for model uncertainty, the 
uncertainty on as is reduced by 30% and the preferred 
value is significantly reduced to as = 0.746 ± 0.017 for 
the homogenous kSZ model. Despite the fact that the 
adopted templa te is the lowest of th e tSZ models shown 
in Figure [8] the Sehgal et al. (20091 template, the value 
of as inferred by the SPT data using this model is lower 
than that favored by WMAP. Additional SPT data will 
soon determine if the apparent tension between as in- 
ferred from SZ and primordial CMB measurements is ro- 
bust. In any case, improving our theoretical understand- 
ing of both the kSZ and tSZ power spectra is essential for 
fully realizing the potential of SZ power spectrum mea- 
surements to constrain cosmological parameters such as 
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7. CONCLUSIONS 



We have presented the first CMB temperature 
anisotropy power spectrum results from the SPT ex- 
periment. Table [T contains bandpowers based on si- 
multaneous observations at 150 and 220 GHz of a single 
100 deg 2 field for three months during the austral winter 
of 2008. A total of 286k detector-hours of data was col- 
lected during these observations, resulting in a final map 
depth of 18 /iK-arcmin at 150 GHz and 40 /iK-arcmin at 
220 GHz. The SPT 150 GHz temperature maps are cal- 
ibrated to an accuracy of 3.6% through a direct com- 
parison of CMB temperature anisotropy as observed by 
SPT and WMAP5. This absolute temperature calibra- 
tion is extended to 220 GHz with a 7.2% uncertainty by 
an internal comparison of the primary CMB anisotropy 
measured by SPT at 150 and 220 GHz. 

The SPT bandpowers include angular multipoles from 
£ = 2000 to 9500. This broad angular range covers the 
transition from where the primary CMB anisotropy dom- 
inates the measured power to where secondary CMB 
anisotropics and foregrounds become dominant. The 
primary anisotropy observed by SPT is consistent with 
ACDM cosmological models based on lower- CMB data. 
At I > 3000, the SPT bandpowers are used to detect 
an excess of power above the primary anisotropies at 
~50 a at both 150 and 220 GHz after masking bright 
(> 6.4 mJy) sources. 

Most of this power can be explained by a population of 
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Fig. 10. — Comparison of the tSZ power spectrum (at 153 GHz) as predicted by numerical simulations and halo model calculations. 
Note that all curves hav e been re-norm alized to the fiducial cosmology. The thick, black, solid line shows the base template, obtained 
from maps generated bylSehgal et al.|d2009t. The black, d ot-da s hed l ine sh ows the modeled spectrum obtained from maps generated 
by applying the semi-a nalytic model for mtracluster gas of |Bode et al,| |2007) to the halos identified in an N-body lightcone simulation 
(as described in IShaw et al.|2009|. The red, solid line shows the results from maps constructed from an adiabatic simulation produced 
using the Eulerian hydrodynamics! code CART ||Kravtsov et al.||2005|l , while the red, dot-dashed line shows results from maps made 
from the MareNostrum simulations | |Zahn et al.|2UIu| |. The blue, solid line shows the predictions of the |Komatsu &~Se ljak (20021 halo 
model calculation. The histogram shows tne Bl-^1 sensitivity to the thermal SZ signal in each £-band (with arbitrary normalization) The 
grey band illustrates the 68% confidence interval of our theoretical uncertainty (see Figure 



faint dusty, star- forming galaxies (DSFGs), the proper- 
ties of which are discussed in H09. We expect the DSFG 
population to contain a clustered component, which will 
produce a power term of comparable size to the tSZ effect 
at 150 GHz. This clustered component would bias the 
measurement of power from secondary CMB anisotropics 
in single-frequency maps. However, we use SPT's multi- 
frequency data to construct the DSFG-subtracted band- 
powers shown in Figure [5] This frequency differencing 
reduces the clustered point source power to negligible 
levels and significantly reduces the Poisson power. 

Secondary CMB anisotropics, consisting of a linear 
combination of the kSZ and tSZ effects, are detected at 
2.6 a in an analysis of the DSFG-subtracted bandpowers. 
Assuming there is no kSZ contribution, this equates to 
4.2 ± 1.5 /zK 2 of tSZ power at i = 3000. This is substan- 
tially lower than the 8.5 fiK 2 predicted by the homoge- 
neous kSZ and fiducial tSZ models using WMAP5 cosmo- 
logical parameters. The lower than expected value of the 
measured SZ power may point to either lower values of erg 
or an over-estimate of the modeled SZ power spectrum. 
Adopting a conservative prior on the uncertainty in the 



model predictions for the SZ amplitude, we combine the 
measurement of the tSZ power spectrum with CMB data 
on larger angular scales to improve the constraint on 
the CMB-derived matter power spectrum normalization, 
(78- This results in a decrease in the preferred value to 
<T8 = 0.773 ± 0.025. A major source of uncertainty in 
this analysis is due to the current theoretical uncertainty 
on the amplitudes and shapes of the thermal and kinetic 
SZ effect power spectra. In the absence of this theoreti- 
cal uncertainty, the constraint tightens about a still lower 
value of as — 0.746±0.017. The tension between the pre- 
dicted and measured SZ amplitudes decreases for models 
that predict less kSZ or tSZ power, for instance models 
with shorter reionization epochs or higher energy feed- 
back in clusters. However, the tension remains at some 
level for all the models we considered. It is our hope that 
the progress in observational efforts represented by this 
work will lead to further work in improving theoretical 
models for secondary CMB anisotropies. 

The SPT power spectrum results presented here use 
only a small fraction of the complete data set. The SPT 
is continuing to take data and has already observed 800 
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Fig. 11. — The ID marginalized erg constraints with and without including the SPT DSFG-subtracted bandpowers for three kSZ cases. 
The black lines denote the erg constraints without SPT, while the red lines include SPT's bandpowers. Constraints with the patchy kSZ 
template are shown with a solid line. The results when including only the homogeneous kSZ model are shown with the dashed lines, and 
the results for no kSZ effect are shown with dotted lines. The SPT data tightens the erg constraint in all three cases. 



deg 2 to similar depth at 150 and 220 GHz to the 100 deg 2 
used in this work, and 600 deg 2 of this area also has 
95 GHz coverage. The complete SPT survey is expected 
to cover over 2000 deg 2 with these three frequencies. The 
complete data set promises a > 10 a detection of the tSZ 
power spectrum, and the addition of the third (95 GHz) 
frequency band will allow us to separate the kSZ effect 
from DSFGs and the tSZ effect. Measuring the kSZ 
power spectrum will open a new window onto the reion- 
ization history of the universe. The complete SPT data 
set will also produce a catalog of hundreds of SZ-selected 
clusters, high signal-to-noise images of known clusters, a 
direct detection of gravitational lensing of the CMB, and 
the widest study to-date on the resolved and unresolved 
source populations at these frequencies. 
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APPENDIX 



BANDPOWER COVARIANCE MATRIX ESTIMATION 



The bandpowcr covariancc matrix includes both signal and noise contributions. The signal covariance is calculated 
from simulations. The noise covariance is estimated from the data. We calculate the variance of the mean power 
spectrum using the variance of cross-spectra between independent real maps. With 300 independent observations of 
the same field, these maps are sufficient to generate an accu rate estimate of the co vari ance matrix. Several details of 



our approach are motivated by the analytical treatments of Tristram et al. (2005) and Polenta ct al. (2005), and we 



first review the analytic estimate of the covariance matrix before discussing the estimator 



2(> 



Analytical Considerations 



Following Tristram et al. (2005), wc represent the expected covariancc between two cross spectra as S: 



^AB.CD 
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Our goal is to express this covariance in terms of the noise and signal in the maps, and to compute the magnitude 
of the diagonal elements, as well as the correlation between bandpowers. As the first step, the central term can be 
rewritten as, 



(DfDg ) (D£ B ) (D° d ) = P bk P Vk , 

= Pbk.Pb'k' 
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To simplify equation A3 Trist ram et al.| ( 2005| ) and Polenta et al. (20051 make the following assumptions 



1. Fluctuations in the map, i.e., from CMB anisotropics, confusion limited point sources and noise, are well described 
by a Gaussian random field. 

2. The beams and filtering applied to the data are isotropic; E^ = G k £? k depends only on |k|. 

3. The instrumental noise is isotropic. 

4. The power spectrum, C kl is smoothly varying with k, and changes little over scales comparable to the width of 
the mode coupling matrix. 

Using the assumption that C k and |_E k | do not vary much over small changes in k, these products become: 



,A~B 
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Here C^ B is shorthand for C k 



Sab, the expected cross spectrum between two unfiltcrcd, perfectly beam- 



corrected — though noisy — maps. The additional term is the n oise bia s that exists in the map auto-spectrum. Assuming 
isotropic beams and filtering, we can combine equations |A3| and |A4| to obtain a relatively simple expression: 
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Combining equations |A3| and |A5| yields: 
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C&°)] (K[W]-\ 3)bl . 

To obtain a simplified expression for the magnitude of the diagonal elements of the covariance matrix, one typically 
assumes the mode-coupling matrix is nearly diagonal: M kk i[W] m W2S kk >. In this approximation (_K"[W] ),,, w 

w^Skk' and thus the covariance of any two cross spectra is: 



w^E b 2 S w and M kk ,[W 2 } 
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where v b is the effective number of independent k-modes in each ^-band. For isotropic filtering, v { 

One subtlety in estimating the covariance is the fact that although the noise in each map is independent, each 
map has the same sky coverage. Hence the signal in all maps is correlated. The correct estimator must take this 
correlation into account. Under the simplifying assumption that all maps are statistically equivalent, the correlation 
between two cross spectra depends only on whether the spectra have a map in common (e.g. the cross spectrum D b 2 
is more strongly correlated to the spectrum D b 3 than D b A ). Comparing the covariance of two cross-spectra taken from 
4 different observations: 
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to the covariance of a pair of cross-spectra with a common map, 

^' BC W^c =1(2^ + ^0,), (A9) 



to the covariance of a pair of cross-spectra with two maps in common, 



The degree of correlation is also £ dependent since it depends on the relative signal vs. noise power in the maps. In 
the high-^ regime, where noise dominates the power in an individual map, all cross spectra are nearly independent. 
Conversely all cross-spectra are nearly completely correlated at \ow-£, where the primary CMB anisotropy overwhelms 
the noise. 

Given the assumption of statistical equivalence, we can then compute the expected variance of the mean spectrum 
based on these variance estimates. This is equal to the correlation between any particular cross-spectrum and the 
mean, 



-bb =-bb Wb « — 2 C b + 4C 6 — — 2 + 2-^4 ) • (All) 



This e stimate of the variation of the mean spectrum agrees with the uncertainty estimates given in Polenta et al 
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It should be noted that the noise and filtering of the SPT data are anisotropic. Equation |A7| can therefore be used 
only as a guideline. In the case of anisotropic filtering or anisotropic beams, the number of independent modes per 
bin will be typically smaller than t^ so , since anisotropic filtering will weight different k-space modes unevenly. For 
example, the k-space mask completely eliminates all modes with k x < 1200. The variance in each £-bin increases with 
fewer independent modes. However, even if we account for the effective number of independent modes in an £-bin, 



equation A~7 does not account for the anisotropic nature of the atmospheric noise contribution. 



The Empirical Covariance Estimator 

The existing analytic treatments are not directly applicable to the SPT data due to anisotropies in the noise and 
filters. By the nature of SPT's scan strategy, atmospheric fluctuations preferentially contaminate low k x modes. 
Likewise the filters intended to remove these fluctuations preferentially remove low k x modes. Instead, we have 
designed an empirical estimator which reproduces the analytical results when applied to isotropic data, while accurately 
accounting for the increased uncertainty due to the noise and filtering anisotropics in the actual data. 

The noise covariance m atrix e stimate is divided into two parts, a signal contribution obtained from the Monte Carlo 
simulations described in §3.2.3| and a noise contribution obtained from real single-observation maps: 

C w = Cf b ?' s + C d b T- (A12) 

The signal contribution is straightforward to estimate with an approach similar to the MASTER power spectrum error 
estimator. We use the signal only simulations to obtain an empirical estimate of the sample variance: 



C MC, S = A £MC, s A £MC, s (A13) 

Note that here Ax = x — x is defined with respect to the sample mean. Since the simulations include only CMB 
realizations and point sources in the confusion limit, the simulated signals are essentially Gaussian. Therefore we 
expect the usual sample variance contribution: 
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V(mc, s) =^ • (A14) 



As before, v b is the effective number of independent Fourier-modes in each i?-band. 

The noise contribution is computed from the cross spectra of single-observation maps. We use the following estimator 
for the noise contribution: 
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Here /(n Q b s ) is a correction due to the finite number of realizations. In the limit of many observations this function 
asymptotes to unity; we use 300 observations so this term can be ignored. The first term can be identified as the 
sample variance of the cross spectra. The second term accounts for the additional correlations between cross-spectra 
with a common map. 



We can now calculate the expectation value for the noise component of the covariance estimator defined in A15 
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This is combined with the signal, or cosmic variance, component in equation |A13| to get the expectation value of the 
estimator: 
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This agrees with the analytic estimate (equation All I for the variance of the mean 



Multifrequency Cross Covariances 

A multifrequency data set requires an estimate of both the covariance of the each individual set of bandpowers 
(i.e., both the single-frequency bandpowers, and the cross-frequency bandpowers) and the cross-covariance between 
these sets. We naturally expect signal correlations between different sets of bandpower due to the fact that all three 
sets of bandpowers reported here are derived from the same patch of sky. However we also expect noise correlations 
between the 150 GHz x 220 GHz cross spectrum bandpowers and each set of single-frequency bandpowers since the 
noise uncertainty in the cross-spectrum is entirely due to noise in the 150 GHz and 220 GHz data. Thus we compute 
the cross-covariance matrices, C^ l '^\ where i and j denote one of the three sets of bandpowers: 150 GHz, 220 GHz 
and 150 GHz x 220 GHz: 
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Treatment of Off-diagonal Elements 

Given the finite number of simulations and data maps, we expect some statistical uncertainty in the covariance 
estimate, particularly in the off-diagonal elements. Such uncertainty is not unique to the estimation technique described 
here, rather it is expected for any covariance estimate which is computed from a finite number of realizations. We 
expect the covariance estimates to be Wishart distributed, with n f, s = 300 degrees of freedom. A given covariance 
clement, C ?1 has a statistical variance of: 
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For diagonal elements we expect a standard deviation of y/2/n b s = 1 /-y/150 = 8.1%. In addition, there is a statistical 
uncertainty on the apparent correlation between two bins. If we assume that the true correlation between bins is 
small, then the standard deviation of the apparent correlation between two bins is y/l/n b s or 5.7%. For the choice 
of bin-size, the statistical error on the correlation of any two bins is much larger than the expected correlation (i.e. 
the fractional error on the apparent correlation estimates is greater than 100% even for adjacent bins). For bins that 
are widely separated, these false bin-bin correlations may skew model fitting. Therefore we "condition" the published 
covariance matric es in order to reduce this statistical uncertainty on the covariance matrices. 

From equation A5 we see that the shape of the correlation matrix (i.e. the relative size of the on-diagonal to 
off-diagonal covariance elements as a function of bin separation) is determined by the apodization window through 
the quadratic mode-coupling matrix, M[W]kfc. For the I range considered in this work, the off-diagonal elements of 
this matrix depend only on the distance from the diagonal, |fc — k'\. Therefore we condition the estimated covariance 

matrix, Ckk 1 by first computing the corresponding correlation matrix, and then averaging all off-diagonal elements of 
a fixed separation from the diagonal: 
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The bandpowers reported in Tables [T] and [2] are obtained by first computing power spectra and covariance matrices for 
a bin-width of A t = 100 with a total of 80 preliminary bins. This covariance matrix was then conditioned according 
to equation |A20| before averaging the bandpowers and covariance matrix into the final bands. 

Equation |A5| is based on the assumption that the filtering is isotropic. In order to test the validity of this equation 
for the anisotropic filtering, we perform 10000 simple Monte Carlo simulations. In each simulation a white noise 
realization is subjected to a simplified, though similarly anis otro pic, version of the filtering scheme. The variance of 
the resultant spectra is computed and compared to equation ~KE Though the apparent correlations be twee n all bins 
exhibit the expected 1% scatter, the correlations between neighboring bins are consistent with equation A5 



